RENORMALIZATION AND BLOW UP FOR CHARGE ONE EQUIVARIANT 

CRITICAL WAVE MAPS. 



J. KRIEGER, W. SCHLAG, AND D. TATARU 

Abstract. We prove the existence of equivariant finite time blow-up solutions for the wave map problem from 
— > of the form u{t,r) = Q{\{t)r) + 1^(1, r) where u is the polar angle on the sphere, Q{r) = 2 arctan r 
is the ground state harmonic map, A(t) = t~^~", and TZ{t,r) is a radiative error with local energy going to 
zero as t — ► 0. The number ^ > \ can be described arbitrarily. This is accomplished by first " renormalizing" 
the blow-up profile, followed by a perturbative analysis. 



1. Introduction 

We consider Wave Maps U : M^+^ S"^ which are equivariant with co- rotation index 1. Specifically, they 
satisfy U{t^ lox) = uoU{t, x) for lo G 5*0(2, M), where the latter group acts in standard fashion on M^, and the 
action on S'^ is induced from that on via stereographic projection. Wave maps are characterized by being 
critical with respect to the functional 

U^l {daU,d°'U) da, a = 0,1,2 

with Einstein's summation convention being in force, 9" = m^^dp, ttIq^ = {m"^)~^ the Minkowski metric 
on and da the associated volume element. Also, (•, •) refers to the standard inner product on R'^ if we 

use ambient coordinates to describe u, daU etc. Recall that the energy is preserved: 

E{u)= j {DU{-,t),DU{-,t)) dx = const 

If one instead uses spherical coordinates, and lets u stand for the longitudinal angle, and similarly use polar 
coordinates r,d on R^, we describe the Wave Map by {t,r,0) — > {u{t,r),9), where now u{t,r), a scalar 
function, satisfies the equation 

Ur sin(2M) 

(1-1) -Utt+Urr-\ = , „ 

r 2r^ 

The problem at hand is energy critical, meaning that the conserved energy is invariant under the natural 
re-scaling U U{Xt, Ax) (using the original coordinates and meaning of U). By contrast, the analogous wave 
map problem on R"+^, n > 3 is energy-supercritical in the sense that the natural scale-invariant Sobolev space 
is then , and the corresponding norm ||u|| is not expected to be controlled globally-in-time for general 
initial data, which leads to the general belief that in this case, there should not be a good well-posedness 
theoory for general initial data, irrespective of the target. Indeed, singular wave maps stemming from C°°-data 
have been constructed on background R'^+^ with target in ,Shal| . and with origin R""*"^, n > 4 and for 
more general targets in |Sha2| . 

In the critical case, global well-posedness is expected for hyperbolic targets, while singularity development 
is expected for certain positively curved targets, such as S^. More precisely, numerical evidence in jBij . |Li) 
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strongly suggests singularity development for equivariant wave maps of co-rotation index 1 from R^+^ to 
with smooth data, while wave maps from ]R^+^ to H^, and more generally H*^, fc > 2, are expected to 
preserve the regularity^ of the initial data. Further evidence for possible singularity development up to this 
date specifically, in the co-rotation 1 equivariant case has been recently provided in |Co| . We note that a 
fairly satisfactory understanding has been achieved for small-energy wave maps from K^+^ to general targets 
|'l'ao| ■ | 'l'at| ■ |Kri| . as well as for rotationally invariant wave maps and general initial data |Chr-Tah| . |Strl| . 
In particular, it is known that the latter never develop singularities |Strl| , and that for equivariant wave maps 
of co-rotation index 1, regularity breakdown can only occur in an energy concentration scenario jStr2| . For 
equivariant wave maps, it is known that regularity of the initial data is preserved (see previous footnote) 
provided the target satisfies a geodesic convexity condition |Sh-Tahl . 

Our objective in this paper is to rigorously demonstrate regularity breakdown for equivariant wave maps 
u : — > of co-rotation index 1 with certain regular initial data. More precisely, the data {u, ut) 
will be of class H^^^ x for some 5 > 0. It is well-known that such data result in unique local solutions 
of the same regularity until possible breakdown occurs via an energy-concentration scenario. We note that a 
result of Struwe shows that if the solution is indeed C°°-smooth before breakdown^, such a scenario can only 
happen by the bubbling off of a harmonic map |Str2| : specifically, let Q{r) : — > S"^ be an equivariant 
hamonic map, which can be constructed for every co-rotation index fc S Z (for example, for fc = 1 stereographic 
projection will do). We shall identify Q{r) with the longitudinal angle, as above. Then according to |Str2) . 
if an equivariant wave map u of co-rotation index fc = 1, again identified with the longitudinal angle, with 
smooth initial data at some time > breaks down at time T = 0, then energy focuses at the origin, and 
there is a decomposition 

u{t, r) = Q{X{t)r) + s{t, r), Q{r) a co-rotation fc = 1 index equivariant harmonic map 

where there is a sequence of times ti ^ 0, ti < 0, i — 1,2, . . ., with A(ti)|<i| oo, such that the rescaled 
functions u{ti, converge to Q{r) in the strong energy topology. 

This is borne out by our main theorem. We let Q{r) represent the standard harmonic map of co-rotation 
fc = 1, i.e., Q{r) — 2arctanr. Recall that in the equivariant formulation the energy is 



\^r..2 , 2n , sin (u) 



.2' ' 2r2 
The local energy relative to the origin is defined as 

„ . N f r 1 , 9 sin^(M) 



r dr 



r<t 



r dr 



.2' ' 2r2 

It is well-known that for equivariant wave-maps singularities can only develop at the origin and this happens 
at time zero iff 

lim inf £ioc (m) (i) > 

The following theorem is the main result of this paper. Note that we need to " renormalize" the profile Q{rX{t)) 
by means of a large perturbation (denoted u"^ below). We find it convenient to solve backwards in time, with 
blow-up as t — > 0-f . 

Theorem 1.1. Let > ^ be arbitrary and > 6e sufficiently small. Define X{t) — t^^^^ and fix a large 
integer N . Then there exists a function'^ u'^ satisfying 

u" £ C+^/^-iitn >t>0, \x\< t}), £i^^{u''){t) < (tA(i))-^ | logip as t ^ 

and a blow-up solution u to in [0, to] which has the form 

u{r, t) = Q{X{t)r) + u^r, t) + e{r, t), 0<r<t 

^By this we mean that if initial data have regularity H^'^^ , 5 > 0, the Wave Map can be uniquely globally extended in this 
class. 

^This result most likely can be adapted to solutions of lesser smoothness 

■^We refer to this as an "elliptic profile modifier"; see Section|^for a detailed explanation of this notion. Also, for noninteger 

(3 means C^I3]^f-m 
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where e decays at t — 0; more precisely, 

eet^H^+'-iR^), Stet^'-^H.^-iR'), S,^^{e){t) <t^ as t ^ 

with spatial norms that are uniformly controlled as t — > 0. Also, u{0,t) = for all < t < tg. The solution 
u{r, t) extends as an H^^^^ solution to all of R? and the energy of u concentrates in the cuspidal region 
< J; leading to blow-up at r = t = 0. 

We remark that a somewhat sm'prising feature of our theorem is that the blow-up rate is prescribed. This 
is in stark contrast to the usual modulation theoretic approach where the rate function is used to achieve 
orthogonality to unstable modes of the linearized problem. Heuristically speaking, there are two types of 
instabilities which typically arise in linearized problems: those due to symmetries of the nonlinear equation 
(typically leading to algebraic growth of the linear evolution) and those that produce exponential growth in 
the linear flow (due to some kind of discrete spectrum). For example, the latter arises in the recent work 
on "center-stable manifolds", see |Sch| . |KrSchl| . jKrSch2| whereas for the former see |KrSch3| . Both types 
can lead to blow up. Here we do not have any discrete spectrum in the linearized equation, but rather a 
zero-energy resonance which is due to the scaling symmetry. It is unclear at this point which role (other 
than a technical one) the resonance plays in the formation of the blow-up. Indeed, our approach is really 
non-perturbative as the crucial elliptic profile modifier produces a large perturbation of the basic profile Q. 
The perturbative component of our proof deals with the removal of errors produced by the elliptic profile 
modifier (it is crucial that these errors decay rapidly in time). 

A recent preprint* by I. Rodnianski and J. Sterbenz |Ro-Stj details the construction of generic sets of initial 

data (including smooth data) resulting in blow-up with a rate \{t) ~ for equivariant wave maps from 

to S*^ with co-rotation index k > A. These data can be chosen arbitrarily close to the corresponding 
co-rotation k harmonic map with respect to a suitable norm stronger than || .||_ffi • This latter behavior appears 
specific for sufficiently large co-rotation indices but, due to numerical experiments (e.g., |Bi| ) is not expected 
for the case of co-rotation index one. More precisely, numerical experiments suggest that perturbations of 
the "ground state" Q{r) = 2arctan(r) need to be of a certain size (depending on their profile) to result in 
blow-up. Indeed, our theorem, which is partly based on perturbative techniques, has a non-perturbative fiavor 
in that the elliptic profile modifier cannot be made small at fixed time t = to > and with a fixed profile X{t). 
On a technical level, we remark that the corresponding linearlized operator has zero energy as an eigenvalue 
for fc > 1 but for k = 1 zero energy becomes a resonance (indeed, d\Qk{Xr)\x=i G L^(0, cx)) iff fc > 1 where 
Qk{r) = 2arctan(r'')). 

Our argument is correspondingly divided into two parts: first, we use a direct method, exploiting the 
algebraic fine structure of the system, to find an approximate solution Q{X{t)r) -\-u'^{t,r). Roughly speaking, 
one may think of Ue(., .) as being obtained by a finite sequence of approximations which alternately improve the 
accuracy near the light cone and near the origin. To model the solution near the light cone, one introduces the 
coordinates (a,i) where a = j and reduces to solving an elliptic problem in a by neglecting time derivatives. 
More precisely, one treats time derivatives as error source terms, which get decimated by iterating the elliptic 
construction. Similarly, one improves accuracy near the origin r = by working with the coordinates (i?, t) 
where R = X(t)r, again reducing to an elliptic problem by neglecting time derivatives. This process does not 
lead to an actual solution, as one "keeps losing time derivatives" , which leads to worse and worse implicit 
constants. Thus in a second stage, we construct a parametrix for the wave equation which is obtained by 
passing to coordinates {R,t) where R = X{t)r, r = ■^i"''. This in turn relies on a careful analysis of the 
spectral and scattering theory of the Schrodinger operator which arises by linearizing around Q{r). The 
remaining error is then iterated away by continued application of the wave parametrix. 

A surprising feature of Theorem ll.il is the fact that blow-up may be arbitrarily slow, as we can prescribe 
V arbitrarily large. However, the data leading to this blow-up are not generic, and indeed rather difficult 
to describe. We observe that in solutions blowing up with X{t) ^ t~^-^ were observed numerically, 
corresponding to a transient regime dividing blow-up from global smoothness and scattering. In particular, 
this blow-up rate appears to correspond to a set of initial data of co-dimension one. Our initial data sets 
seem to lie on manifolds of very large co-dimension, which increases with v (we plan to return to a rigorous 
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4 J. KRIEGER, W. SCHLAG, AND D. TATARU 

treatment of this conditional stability issue in a later paper). In particular, it appears unlikely that such 
solutions corresponding to >> 1 would be detected numerically. 

Finally, we note that our technique quite likely lends itself to constructing blow-up solutions for higher 
homotopy indices, too, as well as to problems of a similar nature, such as the critical Yang-Mills equation. 



2. Approximate solutions 

2.1. The elliptic profile modifier. In this section we show how to construct an arbitrarily good approximate 
solution to the wave map equation as a perturbation of a time-dependent harmonic map profile 

Mo = QiR), R = r\{t) 

with the polynomial timescale 

\[t) = r^-'' 

To describe the approximate solution we use the time variable, the variable R which corresponds to the 
harmonic map scale, and the self-similar variable a = r/t which is useful in analyzing the behavior near the 
cone. The only trade off in this construction is that we need to allow singularities of the form 

as we approach the cone. Thus the larger the parameter u, the better the regularity of the approximate 
solutions. 

Theorem 2.1. Let fc G N. There exists an approximate solution U2k for (II. 1() of the form 

u.._i(r,t) = Qmr) + (^i?log(l + R') + O ^Ell^^S^Lt^^ 
so that the corresponding error has size 

'fl (log(2 + fl))^'=- i' 

i2(tA)2fe 

Here the O(-) terms are uniform in < r < t and < t < tg where to is a fixed small constant. 



e2k-i = O 



Remark 2.2. In the proof wc obtain U2fc-i and 62^-1 which are analytic inside the cone and (72+"", respec- 
tively C~2+''~ on the cone, with a good asymptotic expansion both on the R scale and near the cone. 
More precisely, using our notations defined below we have 

U2k-i e Q{X{t)r) + -l-IS^{RlogR, Qk) 

{txy 

while the error satisfies 

*'^2fe-i e ^-l^/^i(i?(iogi?)2fc-i,Q^_,) 

Proof. We iteratively construct a sequence Uk of better approximate solutions by adding corrections Wfc, 

Uk =Vk+ Ufc-l 

The error at step k is 

/ r^9 .^9 1 X sm(2uk) 

ek = {~d^ + dl + -dr)uk 

r Ir'^ 

To construct the increments Vk we first make a heuristic analysis. If u were an exact solution, then the 
difference 

e = u- Uk^i 

would solve the equation 



X / .^9 .^9 1 ^ X cos(2ufe„i) . , sin(2ufe_i) ,^ ,„ , , 

(2.1) i-df + dl + -dr)e sin(2e) + —\^^^(\ - cos(2£)) = Ck-i 
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In a first approximation we linearize this equation around e = and substitute Uk-i by uq. Then we obtain 
the linear approximate equation 



(2.2) [-d? + d^+'-d,.-'^^je^e,^. 

For r <^ t we expect the time derivative to play a lesser role so we neglect it and we are left with an elliptic 
equation with respect to the variable r, 

(2.3) fa.^ + ia.-H^V«e.„„ r«t 



For r t we can approximate cos(2uo) by 1 and rewrite (|2.2() in the form 

-Or - 

Here the time and spatial derivatives have the same strength. However, we can identify another principal 
variable, namely a = r/t and think of e as a function of {t,a). As it turns out, neglecting a "higher order" 
part of Cfe-i which can be directly included in e/c, we are able to use scaling and the exact structure of the 
principal part of ek-i to reduce the above equation to a Sturm-Liouville problem in a which becomes singular 
at a = 1. 

The above heuristics lead us to a two step iterative construction of the VkS. The two steps successively 
improve the error in the two regions r <^ t, respectively r fa t. To be precise, we define Vk by 



1 „ cos 



(2.4) ( + ^dr - ^-^^] = 



respectively 
(2.5) 



-dt +d^ + -dr - \] V2k ^ e^fc_i 
r r'^ ) 



both equations having zero Cauchy data^ at r = 0. Here at each stage the error term is split into a principal 
part and a higher order term (to be made precise below), 

e/c = el + e\ 

The successive errors are then computed as 

e2fc = <^2k-\ + N2k{v2k), e2k+l = e2fc - d^V2k+l + A^2fc+l(w2fc+l) 

where 

/o R^ AT i \ cos(2uo) - cos(2u2fe) sin(2u2fe).-, ^. , cos(2u2fe) . 

(2.6) -A^2fe+i(w) = V + — (1 - cos(2u)) + — {2v - sm(2w)) 

respectively 

, ^ 1 - cos(2u2fc-i) , sin(2u2fe-i) , cos(2u2fc--i) , . , 

(2.7) -N2k{v) = -2 V + (1 - cos(2w)) + (2v - sm(2'u)) 

To formalize this scheme we need to introduce suitable function spaces in the cone 

Co = {(r,t) : < r < i,0 < t < io} 

for the successive corrections and errors. We first consider the a dependence. For the corrections Vk we use 

Definition 2.3. For i g N we let = i if v is irrational, respectively j{i) — 2i^ if v is rational, 
a) Q is the algebra of continuous functions q : [0, 1] — > R with the following properties: 

(i) q is analytic in [0, 1) with an even expansion at 

(ii) Near a — \ we have an absolutely convergent expansion of the form 

oo / j(2i-l) j(2i) 

q = qoia) + J] (1 - a)^'-'^''+^ ^ g2.-i,, (a)(log(l - a)y + (1 - a)2'-+i 92M(a)(log(l - a)y 

4=1 \ 3=0 j=0 



^The coefficients are singular at r = 0, therefore this has to be given a suitable interpretation 
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with analytic coefficients go, Qij 

b) Qm is the algebra which is defined similarly, with the additional requirement that 

qij{l) =0 if i> 2m + 1, odd 

For the errors we introduce 

Definition 2.4. a) With as above, Q' is the space of continuous functions q : [0,1] ^ with the following 
properties: 

(i) q is analytic in [0, 1) with an even expansion at 

(a) Near a = 1 we have a convergent expansion of the form 

oo / i(2i-i) i(2i) 

q = qoia) + ^ " a)^''"'^''-* E fe-i,,(a)(log(l - a)y + (1 - a)^- ^ fej(a)(log(l - a)y 

i=l \ j=0 j=0 

with analytic coefficients qo, qij 

b) is the space which is defined similarly, with the additional requirement that 

= * ^ 2m + 1, odd 

Next we define tlie class of fimctions of R: 

Definition 2.5. S'"'(i?'^(logi?)^) is the class of analytic functions v : [0,00) U. with the following properties: 
(i) V vanishes of order m at R = 
(a) V has a convergent expansion near R = 00, 

v= E c^jR^-^'ilogRy 

a<j<e+i 

We also introduce another auxiliary variable, 

^^•^^ ^- — (Up — 

Since we seek solutions inside the cone we can restrict 6 to a small interval [0, bo] ■ We combine these three 
components in order to obtain the full function class which we need: 

Definition 2.6. a) S "^{R'' {log RY,Qn) is the class of analytic functions v : [0,oo) x [0,1] x [0, 60] —^'M.so 
that 

(i) V is analytic as a function of R,b, 

u:[0,oo)x[0,6o]^Q„ 

(a) V vanishes of order m. ai R = 

(Hi) V has a convergent expansion at R = 00, 

v{R,;b)= E Cij{;b)R''-^\l0gRy 

o<j<e+i 

where the coefficients Cij : [0, bo] — > are analytic with respect to b 

b) I S™'{R'^ {log RY , Qn) is the class of analytic functions w on the cone Co which can be represented as 

w{r,t)=v{R,a,b), ^; e 5'"(i^'=(logi^)^ Q„) 

We note that the representation of functions on the cone as in part (b) is in general not unique since R, a, b 
are dependent variables. Later we shall exploit this fact and switch from one representation to another as 
needed. We shall prove by induction that the successive corrections Vk and the corresponding error terms Sk 
can be chosen with the following properties: For each A; > 1, 

(2.9) V2k-i e -^IS''{R{logRf''-\ Qk-i) 
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(2.10) fe2k-i e ^±^IS\R{\ogR)'^-\ Q'^_,) 

(2.11) e jj-^ISHR^logRf^-', Qk) 

(2.12) t\2k G -^^[lS\R-\\ogRf\Qk) + bIS\R{\ogRf''-\Q'^)\ 
Moreover, for A: = 0, 

(2.13) t^eo^S^iR-^) 

Step 0: We begin the analysis at k ~ 0, where we explicitly compute cq. 
We have 

Co = -uott 

= -\\'{t)\\^Q"iR)^X"{t)rQ'iR) 

R^Q"(R)^^RQ'(R) 



i2 V ^ (l + i?2)2 ^ 'l + i?2 

1 / , 4i? , 2i? 

-(^^ + 1)7^^55^ + ^^ + !)- 



i2 y ^ > (l + i?2)2 ^ 'l + i?2 



With our notations 



t^eo G ^^(i?-!) 

as claimed. It remains to complete the induction step. Hence we assume we know the above relations hold 
up to fc — 1 with A: > 1, and construct V2k-i, respectively V2k, so that they hold for the index k. 

Step 1: Begin with e2k-2 satisfying H2.12|l or H2.13|l and choose U2fc-i so that (|2.9(l holds. 
If = 1, then define eg := cq. If fc > 1, we define the principal part e§^_2 of e2fe-2 by setting 6 = 0, i.e., 

e§fc_2(i?,a) := e2k~-2{R,a,0) 
For the difference we can pull out a factor of b and conclude that 

t'el,_2&j^^^[lS\R'H^ogRr-',Qk-^)+IS\Ri\ogRr-^Q',^,)] 

C j^IS\Ri\ogR)'''-\Q',_,) 

which can be included in 62^-1, cf. (|2.10(l . 

We define U2fe-i as in 12.4|l neglecting the a dependence of e2j,_2. other words, a is treated as a parameter. 
Changing variables to R in (|2.4(l we need to solve the equation 

{t\)'Lv2k-l = t\l^_2 G ^-l^/5l(i?-l(logi?)2'=-2, g,_i) 

where the operator L is given by 

2 , 1 o cos(2z.o) _ .2 , 1 o 1 l-QR' + 

^-^R+R^^-^^-^R+R^^-J^ (l+i?2)2 

Then H2.9|) is a consequence of the following ODE lemma. 
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Lemma 2.7. Let k > 1. Then the solution v to the equation 

Lv = fe S\R-\\ogR)^''~^), v{Q) = v'{0) = 

has the regularity 

V e S^{R{\ogRf''~'^) 

Proof. Since / is analytic at with a hnear leading term, one can easily write down a Taylor series for u at 
with a cubic leading term. 

It remains to determine the asymptotic behavior of v at infinity. For this it is convenient to remove the 
first order derivative in L (to achieve constancy of the Wronskian). Thus, we seek a solution of 

Ly/Rv:^y/Rf, L ^ dl ^ ' ^ 



4i?2 ^ (1 + i?2)2 
We use this fundamental system of solutions for L: 

^(m-^L_ .,^. -l + ARHogR + R^ 

Their Wronskian is W{(j), 6) = 2. This allows us to obtain an integral representation for v using the variation 
of parameters formula, which gives 

V = -R-h{R) / (l){R')VR'f{R') dR' - / 0{R')y/R'f{R') dR' 

Carrying out the integration shows that the right-hand side grows like R{logR)^''^^ as claimed. □ 
As a special case of the above computation we also note the representation for vi , 

(2.14) VI = J^V{R), V e S^RlogR) 

Step 2: Show that if V2k-i is chosen as above then ()2.1U|) holds. 
Thinking of V2k-i as a function of t, R and a we can write 62^-1 in the form 

e2fc-l = N2k-liv2k-l) + E*V2k-l + E'^V2k-l 

Here N2k-i{v2k-i) accounts for the contribution from the nonlinearity and is given by (|2.6fl . E^V2k-i contains 
the terms in 

dfV2k-iit,R, a) 

where no derivative applies to the variable a, while E°'V2k-i contains the terms in 

i-d^ + d^, + -dr)v2k-i{t,R, a) 
r 

where at least one derivative applies to the variable a. We begin with the terms in N2k-i- We first note that, 
by summing the Vj over 1 < j < 2fc — 2, 

(2.15) U2k~2 - wo G -^IS\R\ogR, Qk-i) 

To switch to trigonometric functions we need 
Lemma 2.8. Let 

ve-^ISHR\ogR,Qk^i) 

{txy 

Then 

sin^^ e -^LS^iRlogR, Qk-i), cosv e /^•'(l, Qk-i) 

[txy 
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Proof. We write 

sinv = vg{v'^) 

with g an entire function. Then it suffices to show that g{v^) G IS'^{1, Qk-i)- We begin with 

e j^IS'iR^ilogR)^, Qk-i) C -^J^IS\{\ogR)^ Qfe_i) 

But a? — R^{t\)~^ and the remaining (iA)^^ together with up to two logi? factors combines to give one h 
factor. We conclude that 

For w G 5*^(1, Qk-i) we evaluate g{w). Since g is analytic we conclude that giw) is analytic in R,b when 
interpreted as 

g{w) : [0,(X3) X [0, 60] ^ Qfe-i 
We consider the asymptotic expansion at i? = 00. Since we have an absolutely convergent asymptotic 
expansion for w and a convergent Taylor series for <? at 0, we obtain an absolutely convergent asymptotic 
expansion for g{w). This gives 

5He50(l,Qfe_i) 

and concludes the proof of the lemma. □ 
Using Lemma [2.81 and 12.15|l we compute 
cos(2uo) - cos(2u2fe-2) = 2cos(2uo) sin^(u2fe-2 - uo) + 2 sin(2uo) sin(u2/c-2 - uq) cos{u2k-2 - uq) 

c -^IS'iR\\ogR)', Qk-i) + j^IS^logR, Qk-i) 



Hence 



2 cos(2wo) - COs(2M2fe-2) 
t ^ V2fc-1 



e ^ (^-^ISH^ogR,Qk-i) + ^IS'{R'{\ogRr,Qk-i)^ .^IS^RilogRf'-^Qk-i) 
C (^/5^(i^-l(logi^)2^ + ^IS^(R(logRf''+\ Q,_i) 

C j^IS^RilogRf^-'Mk-i) 

where at the last step we have pulled a b factor out of the second term. Similarly we have 
2Sni(2u2fe^ 

t (1 - COs(2w2fc-l) 

e ^ (^IS\R-\ Qk-i) + -^IS'iRlogR, Q,_i)) (y-^IS\R{\ogRf^-\ Qk-i[ 
^ ^r^IS''{R-\logR)^''-^,Qk-i) + -^IS'{R{logR)^''-\Qk-i) 



C j^ISHRilogRf^-^Qk-i) 
where we have used a power of b'' to pass to the final inclusion. Finally, 



(2^2.-1 -sin(2z;2,-i)) G ^-^IS%1, Qk-i) (^j^ISHmog Rr'K Qk-i) 

1 



c J^X^IS'iRilogRf'^-^Qk-i) 
C j^IS'iRilogRf^-^Qk-i) 
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This concludes the analysis of N2k-i{v2k-i)- We continue with the terms in E*V2k~i, where we can neglect 
the a dependence. Therefore, it suffices to compute 

Finally, we consider the terms in E°'V2k~i- For 

V2k-iir, t) = j^w{R, a), w e S^{R{\ogR)^''-\ Qk-i) 

we have 

t^E°'V2k-i = JJ-^^ \2kvaWa{R, a) - {v + l)RawRa{R, a) + 2R^^a^^WRa{R, a) + a^^Wa{R, a) 

+ (1 - a'^)waa{R,a) - aWa{R,a)] 
Since Qk-i are even in a we conclude that 

ada, a~^da, (1 - a^)dl : Qk-i ^ Q^i 
Also the R^^ factor simply removes one order of vanishing at R = 0. Hence we easily obtain 

t'E'^v2k-i e ^^IS\R{\ogR)'''-\Q',_,) 

This concludes the proof of ()2.1Uf) . We remark that for the special case of fc = 1, i.e., with vi as in ()2.14|) . 
these arguments yield 

(2.16) t^ei e (U) -2 5-3 (7? log i?) 

Step 3: Define V2k so that (I2.11() holds. 

We begin the analysis with e2k-i replaced by its main asymptotic component /2fe-i at R = oo for 6 = 0. This 
has the form 

2k-l 

t'hk-i = TTTTlfc E 'Z.(«)(logi?)^ e Q'k-i 

which we rewrite as 

^ 2k-l 

t^f2k-i = T^yvik^ E aq3{a){\ogRy 

We remark that H2.16|l implies that fi{a) — (tA)~^alogi?. Consider the equation H2.5|l with /2fe-i on the 
right-hand side. 



[ -df + d', + -dr - ^ j W2k - t'/2fc-l 

rould look for a solution t 
4rT E W^kmogRY 



Homogeneity considerations suggest that we should look for a solution W2k which has the form 

2k-l 

i 

W2k = 



The one-dimensional equations for Wi^^ are obtained by matching the powers of log R. This gives the system 
of equations 



T 



where 



^2 = (j + 1) [((^ + 1)^^(2^ - 1) + «"')^2*'' + - (1 + ^»daWit\ 

+ {3 + 2){j + l){{^+lf+a-^)Wi+^ 
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Here we make the convention that T4^|j, = for j > 2k. Then we solve the equations in this system successively 
for decreasing values of j from 2& — 1 to 0. 
Conjugating out the power of t we get 

{-{dt + ^^^y^) ' + 9'r + ldr-^) Wi,{a) = aq, (a) - F, (a) 
which we rewrite as an equation in the a variable, 

(2.18) L^2k-i).Wi = aqj{a) - Fj{a) 
where the one-parameter family of operators Lp is defined by 

(2.19) = (1 - a^)dl + {a-^ + 2a/3 - 2a)a„ + {-(3^ + /? - a'^) 

We claim that solving this system with Cauchy data at a = yields solutions which satisiy 

(2.20) W4 e a^Qk, j = 0,2k- 1 

To prove this we need the following 

Lemma 2.9. Let < m{j) < j^. Let f be an analytic function in [0, 1) with an odd expansion at and an 
absolutely convergent expansion near a = 1 of the form 



(2.21) 



/(a) = /o(a) + J2 



m(2j-l) 

(l_„)(2.-i)-i ^ /2,_i,„(a)[log(l-a)]-+ 

m=0 

m(2j) 

+ (1 - af^-" ^ /2,,„(a) [log(l - a)Y 



rn=0 

with fij analytic near a= 1. Then there is a unique solution w to the equation 

(2.22) i(2fc-i). w = f, w{0) = 0, daw{0) = 

with the following properties: 

(i) w is an analytic function in [0, 1) 

(a) w is cubic at and has an odd expansion at 

(Hi) w has an absolutely convergent expansion near a = 1 of the form 



(2.23) 



w{a) = wo{a) + 



^(2j-l) 

(l_^)(2.-l).+ i ^ ^i;2,_i,,(a)[l0g(l-a)]^ + 

em 

+ (1 - af^-'+' J2 W2jAa)[log{l - a)Y 

£=0 



with Wij analytic near a = 1 and £{i) = m{i) with one exception, namely £{2k — 1) = m{2k — 1) + 1. // 
however /2fe-i,m(2fc-i) (1) = 0, then £(2k — 1) = m{2k — 1). In that case also W2*;-i.£(l) = if £ > 0, but not 
necessarily when £ = 0. Finally, if f2i-i,j{l) = for all i > k and all j, then also W2i-i,e{l) = for all i > k 
and all I. 

Proof. Denote /3 = (2fc — l)v. Since k>l and v > \, also /3 > ^- Cearly, Lp is the sum of two pieces: that 
part which is homogeneous of degree —2 in a, viz., 

-^^ = + a~^9a - = a~'^da{ada) - 

and the remainder which is homogeneous of degree 0. The equation L^y = has fundamental solutions a 
and a~^. A standard power-series ansatz then leads to fundamental solutions of Lpy = of the form 

Ma) = a(l + 0{a^)), Ma) = a"'(l + 0{a^)) 
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where the O(-) terms are analytic functions of a^. Since our right-hand side / has size 0{a) at 0, this implies 
that we can use the equation to write a complete Taylor expansion for w near 0. Matching coefficients in 
Lfj w — f with 

oo oo 

/(a) = ^(«) = E 

yields 

2j(2j - 2)7«, = (2j(2j - 1) - (4j - I)/3 + P^)wj-i + 
where we take wi ~ 0. The coefficient of Wj is always nonzero; this allows us to successively compute the 
coefficients Wj. The convergence of the series for w easily follows from the convergence of the series for /. 
It remains to study the solution w near a = 1. The behavior of at 1 is well approximated by 

Ll = 2(1 - a)dl + (2/3 - l)da + {(3^0^-1)^ 2(1 - a f+^da [(1 - aj-^+^da] + (/? - /J^ - 1) 
in the sense that 

(2.24) = + (a - 1) [(1 - a)dl + (2(/3 - 1) - a-')da + [a + l)a-^] =: + (a - 1)L} 
The differential operator 

(2.25) 2{l~af+ida[{l~a)-''+ida] 

annihilates 1 and (1 — a)^~^i . Therefore, we seek a fundamental system for y — of the form 



oo 



(2.26) 01(a) = l + ^/i,(l-a)^ 02(a) = (1 - 0)*^+^ [l + ^ AKl - a) 

e=i i=i 
This leads to the conditions, for £ > 1, 

(2.27) Aii(l-2/3)+/3-/32-l = 0, M£+i(^ + 1)(2^ + 1 - 2/3) + (/3 - /J^ - 1)^, = 

(2.28) fii{2(3 + 3) + (3 - 0^ - 1 = Q, fie+i{£ + 1)(2£ + 3 + 2fi) + {/3 - 0^ - = 

Clearly, (|2.28(l always has a solution whereas (|2.27|l requires (i — ^ ^ Z+; in the latter case, the series in H2.26|l 
define entire functions. If, on the other hand, £q :— (i ^ ^ ^ then 0i is modified to 

00 

(2.29) 01 (a) = 1 + ^ Ai^(l - af + Ci 02(a) log(l - a) 

1=1 

with the unique choice ci = —(2/3 + l)^^(/3 — /3^ — 1)m^o- Here (|2.27|) is unchanged and can be solved for 
/i£ up to £ < ^o; for £ > £0 this equation is then modified by the terms from the 02 series (in particular, for 
£ = £q + 1 the choice of ci assures the validity of the equation, whereas for all £ > £q + 1 we can again solve for 
lig). Finally, observe that the same process also leads to a fundamental system for Lp; indeed, the remainder 
(a — 1)L'p in H2.24|l does not change the coefficients of fie+i or fie+i in H2.27|l and H2.28|l . In conclusion, the 
preceding power series argument leads to a fundamental system of Lpy — 0, which we again denote by 0i (a) 
and 02(a). 

Modulo a linear combination of 0i , 02 it suffices to find one solution to the inhomogeneous equation 
Lfsw = f near a = 1. At this point, it will be convenient to write Lp as a Sturm-Liouville operator. Thus, 
we write 

= qi^{a)da{q2{a)da) + qs{a) 

with, cf. (I^TOjl . 

griq2(a) = l-a2, ig^(a) = a"! + 2a(/3 - 1), 93(0) =-/?'+/?- a"^ 

One checks that for a close to 1 the first two equations admit solutions 

q^ia) = (1 - a)-^+^[l + (1 - a)gi(a)], gi(a) - ^(1 - a)-''-^[l + (1 - a)^^)] 
with qi , 52 analytic near a = 1 . The Wronskian can now be computed as 

q2{a)[c^i{a)c^'M " ^'M<l>2{a)] = /3 + 1/2 
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Thus, a particular solution of the inhomogeneous problem is given by 

(2.30) w{a) = (/3 + 1/2)- Vi (a) 02(a')'Zi(«')/(a') da' + (/3 + 1/2)- V2(a) T </>i(a')'Zi(a')/(a') da' 



where ao < 1 is some number close to 1. For the first integral, note that (j)2{a')qi{a') is an analytic function 
in the neighborhood of a' = 1. Let 7 7^ —1 and m be a positive integer. Iterating the relation 

(2.31) (\l - a')^[log(l - a')r da' = ^[log(l - a)]™(l - a)^+i - ^ f\l - a'r[log{l - a')]""' da' 
J a 7+1 7+1 Jq 

shows that each summand on the right-hand side of (|2.21|) . inserted into the first integral in 12.30|l . makes an 

admissible contribution to w in the sense of (|2.23|) (for this it does not matter whether <j> takes the form H2.26|l 

or H2.29|l 'l. The analysis of the second integral in (|2.30l) is again based on H2.31|l provided j ^ k, since then 

7 7^ —1. li j = k, then we encounter 

'(1 - a'y^[log{l - a')r da' = -(m + l)-i[log(l - a)]™+i + C 

I ao 

which explains why we obtain an extra log-factor when j = k. Clearly, if /2fc-i,jn(2fe-i) (1) = then there is 
no extra log-factor and the lemma is proved. In that case also we write, with 

m(2fc-l) 

f^{l-af-i /2fe-i,m(a)[log(l-a)]™ 
the second integral in H2.3U|I as 



<^2(a) / (t)i{a')qi{a')f{a') da = (j)2{a) / (t>i{a')qi{a')f{a') da ~ 4>2ia) / 4>i{a')qi{a')f{a') da 

J ao ao ^ a 

The first term on the right-hand side here is just a multiple of 02(0)7 whereas the second one possesses the 
extra vanishing at a = 1, as claimed. The final claim of the lemma follows similarly. □ 

Before turning to the proof of H2.20|l in full generality, we first discuss the special case fc = 1 . This will also 
serve to explain how the algebra Qk arises at all in the iteration. If = 1, then (|2.18|) reduces to the system 

Wlia) = a, W^{a) = -{v{v + 1) + 0-^)1^2^(0) - [a'^ - [v + l)a)daWl{a) 

due to t^fi{a) — (iA)-^a log i?. In view of the solution formula l|2.3()(l with (3 — provided v — ^ ^ Z+, 

Wl{a)^g^[a)+gi{a){l-aY+"- 

W^{a) = ho{a) + hiia)il - a)"+3 + /i2(a)(l - a)'^+3 log(l - a) 

where gj{a), hj{a) are analytic around a — \. Note that the term (1 — 0)^+5 log(l — a) appears in W2 due to 
daW^. Similarly, if - i e Z+, then 

W^{a) = 50(a) + gi(a)(l - aY+"^ + 52(a)(l - aY+"^ log(l ~ a) 

2 2 
W!^{a) = ho{a) + {1- ar+i ^ /i,+i(a)[log(l - a)]' + (1 - 0)^''+^ ^ /i,+4(a)[log(l - a)]^ 

(=0 i=0 

with analytic gj, hj. The terms involving the (1 — a)'^'^+^ factor in W2 are due to the modified (f>, see ()2.29|l . 
Thus, we see that in all cases W2 G Qi for j = 0,1 and a near 1. 

We now continue with the proof of H2.2Q() for general k. At first we consider the easier case when v 
is irrational. We apply the lemma in (|2.18|) using for the right-hand side the fact that q2k-i S Q'k-i- This 
implies that the coefficient of (1 — vanishes at a = 1 . The lemma gives a similar expansion 

for W"^^'"^ with the required vanishing conditions. Hence W^^fc ~^ ^ "^ith one extra (1 — 0)'^'^"^)'^+ 2 term 
(this is the W2fe-i,o(l) 7^ statement of the lemma) - we refer to this as the "free term" in what follows. 

Next we reiterate the argument for the remaining W^fc which solve (|2.18|) . At each step we have to compute 
Fj, see (|2.17|l . Since W^^^ and W^^"^ have an odd Taylor expansion at beginning with a cubic term, it 
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follows that Fj has an odd Taylor expansion at beginning with a linear term. The expansion of Fj around 
a = 1 is similar to the one for W^^^ except that one (1 — a) factor is lost in the "free term". For j = 2k — 2 
this produces the term (1 — a)^^'^~^^''+5 log(l — a) in VK^fc ^tc. At the conclusion of the iteration we have 
gained at most 2k — 1 logarithms in the free term for the VFl^.'s. Then (|2.20() follows. 

Next we consider the case when i/ is rational. This is more difficult since now the term (1 — a)^^'^"-'^)''"^ 
can also arise in expressions of the form 

/2,_i,™(a)(l - a)(2^-i)''-^[log(l - a)r or /2,;„.(a)(l - a)2^'^[log(l - a)]™ 

using the notations of the lemma. This leads to more logarithms than in the irrational case. The first term 
above will be of interest if 2{k — j)^ is an integer, while the second needs to be considered if (2fc — 2j — — ^ 
is an integer. The worst case is j = — 1. Then we can have m(2k — 2) logarithms in the second term above, 
while 2k more logarithms are produced by the 2k applications of the lemma. Hence we need the relation 

m{j) > m{j - 1) + j + 1 

which is verified e.g. by rn(j') = (we pick rij — 2j^ because of j = 1, see above). 

We cannot use 'W2k for W2fc due to the singularity of logi? at i? = 0. However, we define instead 



V2k " 



In doing this we add an additional component to the error. This is large near i? = 0, but this is not so 
important since the aim of this correction is to improve the error for large R. Since — R^/{tX)'^, pulling a 
cubic factor out of the VF's it is easy to see that H2.11|l holds. 

Step 4: For V2k defined as above we show that the corresponding error e2k satisfies H2.12II . We can write 
e2k in the form 

t'esfe = {e2k-i - e°k-i) + (^e^fc_i - (^-d? + + ^dr - V2k^ + t^N2k{v2k) 
where N2k is defined by l|2.7|l and 

„ 2/c-l 

^ ^ j=o 

We begin with the first term in e2k, which has the form 

t^{e2k-i - e"2k-i) e [lS\R-\logRf', Q^i) + MS'iRilogRf^-', Q^i)] 

The second term is contained in the second term of (|2.12(l . It remains to show that 

(2.32) /5l(i^-l(logi^)2^ Q',_,) c IS\R"\\ogRf\ Qk-i) + bIS\RilogRf''-\ Q',_,) 

For w e /S'i(i?-i(logi?)2'=, Q'i^_^) we write 



Then 



w = (1 — a^)w + R^w 

itxy 

(1 - a^)w e IS\R-\\ogRf\ Qk^i), j^R^w E bl S\R{log R)^''-\ Q'^_^) 

{txy 



as desired. The second term in e2k would equal if we were to replace i log(l + R^) by log R in both e^/. and 
V2k- Hence the difference is obtained when we replace the derivatives of i log(l + R^) by derivatives of logi? 
in the expression 



-d^ + c^. + U) V2k = [-d^, + dl + ( 5: Wi,ia) (i log(l + R') 
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Computing these differences one finds that the second term in e2k is a sum of expressions of the form 



i / J n 



Since are cubic at it foUows that we can puU out an a factor and see that this part of the error is in 



' 2k 

1 



which is admissible by (|2.32|) . 

FinaUy, we consider the nonhnear terms in A^2fe- Again the a, h dependence is uninteresting since Qj, is an 
algebra. We shall use that 

1 o 

W2/C-1 - Mo e 7TTT? ■IS%R\ogR,Qk) 

[txy 

By Lemma 12.81 for the linear term we therefore have 

2 1 - COs(2u2fc_l) 
^^2 ^2fc 



C (^IS'{R-\\ogR)''^-\ Qk) + j^IS'iRi\ogR)''^-\ Q, 

For the quadratic term we obtain 

2sin(2u2fe-i) , , 

t r-^ (1 - COs(2w2fc)) 

2r^ 

^IS\R-\Q,) + j^IS'{RlogR,Qk)] (— ^ I S%R' {log Rf''-'),Qk) 



i?2 y ^ (tA)2 ' '-""7 V(tA)2fe+2 

C {-^^^ISHR'ilogRr-^ Q,) + ^^±_IS'iR^logRr-\ Q,: 

1 



C {lS\R-\logR)^\ Qk) + bIS\R{\ogRf''-\ Qk)) 



Finally, the cubic term is 



^2COs(^^^2.2.-sin(2.2.)) G ^((^/^^(i?^(logi?)^'=-),Q.: 



1 1 



^ (txr itX)*k+, IS'(R'i^ogRr'^Qk)) 



b 



C ^^/^\i?(logi?)2'=-i,Q',)) 
This concludes the proof of Theorem 12.11 □ 

3. The perturbed equation 

We now need to complement the approximate solution found in the first section to an actual solution. 
The mechanism for achieving this will rely on the construction of an approximate parametrix for a suitable 
wave-type equation. We now set about deriving this equation: we seek an exact solution of the form 

u{t, r) = U2fc_i(i, r) + e(i, r) 
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where U2fc-i is as in the previous section and e will be obtained by means of an iteration procedure. To 
motivate this procedure, note that we need to solve the following equation, see (jOJ, 

/o i\ 1 cos(2(3(Ar)) , , 

(3-1) - ett + err + -£r , ^ = ^2fe-i(e) + ea/c-i 

r 

where A^2fc-i is defined in (|2.6|) but with U2k-2 replaced by U2fc-i- 

In order to remove the time dependence of the potential in Ij^.H . we now introduce new coordinates: first, 
the new time is to satisfy the relation 

d_ _ 1 d 
~ X(t)di 

Specifically, we may put t — — \(s)ds+ ^ ~ h^ '^' Thus, the singularity now corresponds to r = oo. Next, 
introduce the new dependent variable v(t,R) :~ e(t(T), A^^i?), where we now understand A as a function of 
r. Then we have 

—V = t'{T)et{t{T),X-^R) - ^Rer{t{T), \-'R), —v = A-i£,(t(T), A-^i?) 

This entails that 



{dr + ^Rdn)v = X'^Mt), X-^R) 
A 



From here we get 

(dr + ^Rdafv = (dr + ^Rdn)[X-ht] = X'hu - ^et - X'^u - ^d^v - [^^Rdnv 
We conclude that we may recast the wave equation (I3.1() in the following way: 

-[{d, + ^Rdnf + ^{dr + ^Rdn)]v + {dl + ^dn £^![^|^)„ ^ ^[A^2.-i(s) + e.^-iMrl A'^i?) 

In order to turn the above second order elliptic operator in R into a selfadjoint operator we introduce the 
new variable e(r, R) :— R2v{t, R). This leads to 

(dr + ^RdR)e = RHdr + ^RdR)v + lR^^viT,R) 

{dr + ^Rdnfe = Ri (9. + ^i?^^)^^ + + ^Rdn)v + + 

One checks that 



as well as 



R^dr + -fRdRYv + R-^^id, + -fRdR)v = {dr + -fRdnYe^ i^T^ 2^^V 
^t'n^ • cos[2Q(i?)]^ 



4i?2 (l + i?2)2^ 

Combining these observations with (|3.1|) . we now obtain the wave equation 
(3.2) (^-{dr + ^Rdnf + \{^f + \dr{^)^ e - Le = X'^R^ (n2u-i{R-"^ e) + e^k- 

where 

(3-3) C:^-dl + ^- JY^J^ 

Equation (|3.2I) is the main equation which we need to solve in this paper. As a first step, in the following 
section we will carefully analyze the spectral properties of the underlying linear operator C. 
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4. The analysis of the underlying strongly singular Sturm-Liouville operator 



The goal of this section is to develop the scattering theory of £ from We start with the basic^ 

Definition 4.1. Let 

^0--^-^ + ^, C:^Co- j^^, =:Co + V 
be half-line operators on L^(0,cxd). They are self-adjoint with the same domain, namely 

Dom(/:)=Dom(/:o) = {/ei'((0,cx3)) : /, /' £ Aaoc((0, oo)), £0/ e i'((0, oo))} 

It is important to realize that because of the strong singularity of the potential at r = no boundary 
condition is needed there to insure self-adjointness. Technically speaking, this means that Co and C are in 
the limit point case at r = 0, see Gesztesy, Zinchenko jGesZin| . It is worth noting that the potential is 
critical with respect to this property — any number smaller than | leads to an operator which is in the limit 
circle case at r = 0. We remark that Cq and C are in the limit point case at r = 00 by a standard criterion 
(sub-quadratic growth of the potential). 

Lemma 4.2. The spectrum of C is purely absolutely continuous and equals spec(>C) — [0,oo). 
Proof. That C has no negative spectrum follows from 

^3/2 

(4.1) C(t,o=0, Mr)^T-^ 



with (/)o positive (by the Sturm oscillation theorem, see | DS p. And since (po ^ i^((0, 00)), zero is not an 
eigenvalue. The pure absolute continuity of the spectrum of £ is an immediate consequence of the fact that 
the potential of C is integrable at infinity. □ 

Since <po ^ i^((0,cxD)), one refers to zero energy as a resonance. Heuristically speaking, this notion can be 

thought of as follows: by inspection, Cor~2 = Q and Cor^ = 0. A "generic" perturbation £ — Cq V with 
~ ... . ..13 

V bounded, smooth, and nicely decaying, will have zero energy solutions that behave just like r~2 and , 

respectively. However, in some cases V is such that these two £0 solutions will be "in resonance" and produce 

a globally bounded zero energy solution of £ which behaves like r^ around zero and around infinity — 

just like (1)0. 

For the parametrix construction in the following sections the relevance of the zero energy resonance lies 
with the singularity of the spectral measure of £ at zero energy. Indeed, for £0 the density of the spectral 
measure behaves like ^ as ^ ^ 0, whereas for £ we will show that it behaves like (^ log^^)~^ as ^ ^ 0. We 
now briefly summarize the results from |GesZin | relevant for our purposes, see Section 3 in their paper, in 
particular Example 3.10. 

Theorem 4.3. a) For each z £ C there exists a fundamental system <j)[r, z), 9{r, z) for C — z which is analytic 
in z for each r > and has the asymptotic behavior 

1 1 _i 

(4.2) (j){r,z)^r^, 6{r,z)^-r 2 as r ^ 

In particular, their Wronskian is W{9{-, z),(j){-, z)) — 1 for all z £ C We remark that (j){-,z) is the Weyl- 
Titchmarsh solution^ of C — z at r — 0. By convention, (f>{-, z), 9{-, z) are real-valued for z G R. 

b) For each z £C, Im z > Q, let ip^{r, z) denote the Weyl-Titchmarsh solution of £ — z at r = 00 normalized 
so that 

J 1 J 
'0"'"(r, z) ^ z~4 e*^ as r^oo, Imz2>0 



^In this section we use the variable r > for the independent variable. The reader should note that this now plays the role 
of R in the previous section. 

^Our (f>{-, z) is the ^(z, ■) function from IGesZinj where the analyticity is only required in a strip around M - but here there is 
no need for this restriction. 
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If ^ > 0, then the limit ip'^(r,^ + iO) exists point-wise for all r > and we denote it by ijj~^{r,(,). Moreover, 

define ^+(-,^). Then xp^{r,^), ?/;^(r,^) form a fundamental system of C — with asymptotic 

^ 1 

behavior '\li'^{r,S) ^ e^^^^'' as r ^ oo. 

c) The spectral measure of C is absolutely continuous and its density is given by 

(4.3) p(e) = ilm m(e + iO)xK>o] 

with the "generalized Weyl-Titchmarsh" function 

^ ' M/(^+(.,0>(.,0) 

d) The distorted Fourier transform defined as 

: / ^ /(O = lim [\{r,Of{r)dr 
is a unitary operator from L^(R+) to L^(M+,p) and its inverse is given by 



T-':f^ f{r) = lim / 0(r, 0/(0 p(0 

Here lim refers to the L^(M+,p), respectively the L^(M+), limit. 

Needless to say, part b) above has nothing to do with |GesZin| and is standard. Most relevant for our 
computations are H4.4|l (which is formula (3.22) in ,GesZinj), as well as the Fourier inversion theorem in this 
context (see Theorem 3.5 in |GesZinj '). 

Theorem 14 .31 of course also holds for Cq instead of £. In that case we have a Bessel equation with solutions 

(4.5) 0(r; z) = 2z-^^\^^\Ji{z^^\) 

e{r- z) = '^z^/^r^/^[-Yi{z^''^r) + -k'^ log(z) Ji(z1/V)] 

^(r-z) = z'l^r^l^\-Y^{z'l^r)^i,h{z^l^r)\ = z^l^r^l^iH'^l\z^l^r) 
= Q{r\ z) + m{z)(j){r] z) 

to(z) = ^z[i-7r"Mog(z)], zGC\M+ 

The last formula shows that for strongly singular potentials the Weyl-Titchmarsh function ceases to be 
Herglotz, see ' (jCsZinl for further discussion. Although we shall make no use of these formulas for Cq, the 
reader should note the similarities between the asymptotic expansions on 0, 9 and ip'^ we derive below and 
the classical ones for the Bessel functions, cf. jWat| . 



4.1. Asymptotic behavior of (j) and 6. Beginning with two explicit solutions for Cf — 0, namely 

a I \ 1 — 4r^ log r — r'* 
iTT^' "^'^ ^ 2r^(l + r2) 

we shall construct power series expansions for and from (|4.2|) in z G C when r > is fixed. 



<^o(0 = -r-^^ ^o(0 = ^ ^ = r- 2 (1 _ r^)l2 _ 20o(r) logr 



Proposition 4.4. For any z e C the fundamental system cf)(r,z), 0(r,z) from Theorem \4-.!J\ admits absolutely 
convergent asymptotic expansions 

oo 

((>{r, z) = Mr) +r^^ Xl^''^^)^'^^^^^) 

OG 

0{r,z) ^ r-i(l -r^ - '^{r^zy9j{r^)'^ /2 - {2 + z/i)(j){r, z)logr 

3 = 1 
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where the functions (f)j, 9j are holomorphic inU = {Rew>— 5} and satisfy the hounds 

3C^ 1 
\(l)j{u)\ < (j_i)j lQg(l + l"l)' \Mu)\>^iogu 

Mu)\<C\u\, \ej{u)\<-^a + \u\), ueU 



Furthermore, 

(4.6) Mu) 



— jlogM+i + 0(u ^ log u) as « — > 00 
-f + f^+0(u3) as u^O 

Proof. We begin with ^. We formally write 

00 

(l){r,z) =r-i^z^fj{r) 
j=o 

This becomes rigorous once we verify the convergence of the series in any reasonable sense. The functions fj 
should solve 

C{r-ifj)=r-ifj.u /_i = 
The forward fundamental solution for C is 

H{r,s) = ^{Mr)eo{s) - Ms)Oo{r))l[r>s] 

Hence we have the iterative relation 



fj{r) = I I ris-HMr)Oo{s) - Ms)Oo{r))fj-i{s)ds, fair) 



1 _ 1 



Using the expressions for ^0 we rewrite this as 

i- / ^ r r^{-l + 4s'^logs + s^)-s'^{-l + 4r^logr + r'^) . , , , 
^^^'^ = y„ 2,s(l + r^)(l + .^) 

We claim that all functions fj extend to even holomorphic functions in any even simply connected domain not 
containing ±i, vanishing at 0. Indeed, we now suppose that fj-\ has these properties and we shall prove them 
for fj . Clearly, fj extends to a holomorphic function in any even simply connected domain not containing ±?' 
and 0. We first show that at there is at most an isolated singularity. For this we consider a branch of the 
logarithm which is holomorphic in C \ R+ and show that fj{r + iO) = fj{r — iO) for r < 0. Disregarding the 
terms not involving logarithms, we need to show that for any holomorphic function g we have 

(log s — log(r + iO))g{s)ds = / (log s — log(r — iO))g{s)ds 

Jo 

This is obvious since for s < we have 

log(s + iO) — log(r + iO) = log(s — iO) — log(r — iO) 

The singularity at is a removable singularity. Indeed, for s close to we have |/j-i(.s)| ^ |.s| which by a 
crude bound on the denominator in the above integral leads to |/j(r)| < |r| (again with r close to 0). This 
also shows that fj vanishes at 0. 

The fact that fj is even is obvious if wc substitute 2 log s and 2 log r by log respectively log in the 
integral. This is allowed since due to the above discussion we can use any branch of the logarithm. Indeed, 
denoting /j_i(s^) = fj-i{s) the change of variable =v yields the iterative relation 

(4.7) f,{u)=j^ 4.(1 + .)(! + .) 



1 + u 
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Next, we obtain bounds on the functions fj. To avoid the singularity at —1 we restrict ourselves to the region 
U — {Re It > — We claim that the fj satisfy the bound 

\f,{u)\<j—^\u\Hog{l + \u\) 

The kernel above can be estimated by 

u(— 1 + 2i;logt' + v'^) — v{—l + 2wlogit + m^) 



We have 



which yields 



2v{l + u){l + v) 



l/o(")l < 3 



1+ u 



\fi{u)\<3C\u\ 



\u\ 



1+X 



dx = 3C|u|log(l + 



From here on we use induction, noting that 

"I 



x^"^ log(l + x)dx < -\u\^ log(l + |-u|) 
J 

To conclude the proof, we note that the functions are given by 4>j{u) — u^^ fj{u) and satisfy the desired 
pointwise bound. Finally, (|4.6|l follows by an asymptotic evaluation of the explicit integral ()4.7|l with j = 1, 
which we leave to the reader. 

The argument for the function 9 is similar. The ansatz 



r,z) 



OO 

'K^ " " 5Z ^'Sj W) /2 - (2 + ^/4)0(r, z) log 7 



= (l - - ^ z-'5,(r))/2 + (2 + z/4){Mr) - ^'^"'/j W) l^S' 

yields a recurrence relation for the gj via {C — z)9 = 0. Indeed, for j = 1, 

£{r-^gi{r)) = 6*0(0 " /:(^'/>o(f ) logr + 4r"5/^(,^) logr) 



r 2 



r2(3 + r2) 8 



hir) + -mr) 



where the important fact is that the quantity in brackets is even analytic around and vanishes at 0. A 
similar computation yields for j > 2 

The same considerations as in the case of f j show that each gj is an even holomorphic function in any even 
simply connected domain not containing ±j. Also, the same bound for the fundamental solution for C leads 
to |s'i('')| < C'r"' and more generally, for j > 2, 

\9Ar)\<J—Yyr''{l + r') 
The proof of the proposition is concluded. □ 

Remark 4.5. The logarithmic behavior of (t>i{u) for large u is inherited by (/)(r, indeed, suppose that 
1 ^ ^ > and r = (5^^ 2 where (5 > is small. Then the proposition shows that 

'/'('^i C) ^ '^^^ logr 

The size of 5 here only depends on various constants in the expansion of (p and is thus itself an absolute 
constant. We remark that the appearance of the logr term is a specific feature of £ — it does not occur for 
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£o) see (|4.5|l — indicative of the fact that £ is a long range perturbation of Co- We shall see later that the 
logarithm in cj) produces crucial logarithmic factors in the small ^ asymptotics of the spectral density of £, 
see Proposition 14. 71 below. 

We note that although the above series for (j) converges for all r, z, we can only use it to obtain various 
estimates for in the region |z|r^ < 1. On the other hand, in the region ^r^ > 1 where z = ^ > 0, we will 
represent in terms of ip^ and use the ip^ asymptotic expansion, described in what follows. 

4.2. The asymptotic behavior of ip^. The following result provides good asymptotics for in the region 

Proposition 4.6. For any S, > 0, the solution ^p^{-, ^) from Theorem \4-.y\ is of the form 

iP+{r,^)=C^e'''^^crHKr), > 1 

where a admits the asymptotic series approximation 

3=0 

with zero order symbols ipj'{r) that are analytic at infinity, 

sup \{rdr)''ipf{r)\ < oo 

r>0 

in the sense that for all large integers jo, and all indices a, (3, we have 



sup 

r>0 



[rdrTiqd.f a{q,r)-y2q-'yjUr) 



1° 



for all q > 1 . 

Proof. With the notation 

we need to solve the conjugated equation 



2^^rf, ^ ( + — - ) fo 



(4.8) (-d^^ - 2^^^^r + ^ - (YTT^ ) -(^^^^) - 
We look for a formal power series solving this equation, 

oo 

(4.9) Er*/^« 
This yields a recurrence relation for the /j's, 

3 

4^2 ^ ^2 

which is solved by 

Extending this into the complex domain, it is easy to see that the functions fj are holomorphic in C \ [—i, i]. 
They are also holomorphic at oo, and the leading term in the Taylor series at oo is r^^. At 0, on the other 
hand, fj are singular. The worst singularity is of power type, namely ; however, weaker terms contain 
logarithms and powers of logarithms so it is not easy to obtain a complete expansion. Instead we contend 
ourselves with a weaker estimate, namely 

\{rdr)''fj\<Cjkr-^ Vr>0 

which is easy to establish inductively. The functions 
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now satisfy the desired bounds due to the bounds above on fj. 

Unlike in the expansion for small r, here we make no effort to obtain a uniform estimate on the size of the 
derivatives of ipj'. This is because we do not expect the formal series (|4.9|l to converge, on account of the fact 
that derivatives are lost in the iterative construction of the fj 's. Instead we can construct an approximate 
sum, i.e., a function <7ap{q, r) with the property that for each jo > we have 

JO 

(4.10) (r-a.r(ga,)'^[a,p(g,r)-^g-V+M] I^c^MoI^'"-' 

The construction of aap{q,r) is standard in symbol calculus; indeed, we can set 

oo 
3=0 

where 6j — > sufficiently fast and x is a cut-off function which vanishes around zero and is equal to one 
for large arguments. The bound H4.1(}|l implies that aap{r^i,r) is a good approximate solution for H4.8|l at 
infinity, namely the error 

satisfies for all indices a,P,j 

\{rdrnqdg)^e{q,r)\ < c^^p,,r-\-^ 
To conclude the proof it remains to solve the equation for the difference ai = —a + aap, 

-dr - 2^e^a. + A _ ^^^^ <yiHKr) = eK^r) 

with zero Cauchy data at infinity. We claim that the solution ai satisfies 

\{rdrnqdgy'aAq,r)\ < c^,p.,q-^, J > 2 

Note that this finishes the proof by defining a — aap — cri . A change of variable allows us to switch from the 
pair of operators {rdr,qdq) to (rdr,S.d(^) with comparable bounds. We rewrite the above equation as a first 
order system for (fi,U2) = {o'l^rdrai): 



dr 

Then we have 
which gives 
and by Gronwall 



dr 

^\v\ > -C(r->|+r|e|) 
dr 



Hr)\ < I (^) s\eis)\ds 

Then for large j we have 

(4.11) |e| < C^r~^-^ => \v\ < C^r-^ = q~^ 

To estimate derivatives of v we commute them with the operator. For derivatives with respect to r we have 

'^^''^^ ( S )-( ^ - ^ ) ^^'^^ ( S ) = (l: - 4^ l)iZy{ rdXe) 

But the right-hand side is bounded by r^^^^ from the previous step and the hypothesis on e, therefore as 
above rdrV is bounded by r^^ . 
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We argue similarly for the ^ derivatives. We have 



re 



The only difference is in the first term on the right, for which we write ^2 = r and we use the decay 
property of v with j replaced by j ' + 1: 

as desired. Finally, higher order derivatives are estimated by induction using the above arguments at each 
step. □ 

4.3. Structure of the spectral measure of C. We begin by relating the functions 0, 9 and -0^- By 
examining the asymptotics at r = we see that 

(4.12) W{e,(j))^l 
Also by examining the asymptotics as r 00 we obtain 

(4.13) W{il^+ ,^-) = -2i 

Hence we can express the C ~ solutions in cither the (j), 9 basis or the tp^ basis. On the other hand, 6 
are real-valued while the real and imaginary parts of are equally strong. Hence the two bases are quite 
separated. These are the main ingredients of the next result. 

Proposition 4.7. a) We have 

(4. 14) 0(r, = a{i)i;+ (r, + a(0^+(r,0 
where a is smooth, always nonzero, and has size^ 




Moreover, it satisfies the symbol type bounds 

\m)''aio\<ck\am ve>o 

b) The spectral measure p{^)d^ has density 

TT 

and therefore satisfies 

Proof, a) Since (p is real- valued, due to H4.13|l . the relation (|4.14l) above holds with 

We evaluate the Wronskian in the region where both the ilj^{r,£^) and 4>{r,£,) asymptotics are useful, i.e., 
where r^^ « 1. By Proposition 14 . 41 we obtain that both (t>{£,~^,£,) and {rdr4>){^~^ ,£,) can be expressed in the 
form C^/(^~^) with / holomorphic and satisfying 

|/(u)|<log(l-H|«|), Reu>^ 



X 6 means that for some constant C one has C < b < Ca 
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On the other hand, it follows from ProDosition l4 . 61 that both , £,) and (r3r'0''')(C~^ j '^^^ be expressed 

in the form ^~^h{£^~^) with h satisfying symbol type bounds 

\{rdr)''h{r)\<Ck 

Combining the two expressions above, it follows that a is a sum of terms of the form with 
/, h as above. The bounds from above on a and its derivatives follow. 

It remains to prove the bound from below on a, which is more dehcate. By (j^^ we have 

Since (f) is real-valued, this gives 

Im OW{q^{-,0, (-,0)] = -5.0(r, 

which implies that for all r we have 

We use this relation for r = S£^~2 with a small constant S. Then by Proposition 14.41 we have 

\drcf>{r,0\>r-hog{l + r^) 

while by Proposition ^21 

This give the desired bound from below on a. 

b) By (|4.12|l we can express in terms of 6 and cj) by 

V;+ = -<j>W{ip+,0) + 9W{^J+,<j)) 

Since both (f) and 9 are real-valued, by inserting this into (|4.13|) we obtain the relation 

lm{W{^+,e)Wi^p~',(j))) = -1 

Inserting this in the expression for the spectral measure 1)4. 3|l and taking (|4.4|) into account we obtain 

1 lm{W{4'+,0)Wi^-,^)) _ 1 ,+ _ 1 



as desired. □ 

5. The transference identity 

Returning to the radiation part e in l|3.2|l . the idea is to expand it in terms of the generalized Fourier basis^ 
4>{R,^) associated with the operator C — ~dj^ + — ^^^^ap , i.e., write 



Jo 

and deduce a transport equation for the Fourier coefficients x(r, ^). The main difficulty in doing this is caused 
by the operator Rdu which is not diagonal in the Fourier basis. Our strategy for dealing with this is to 
replace it with 2^9^ modulo an error which we treat perturbatively. The operator Rdu — 2^9^ is natural since 

it annihilates the expression e'^^^ arising in the asymptotic expansion of <j){R,^) for large R. Consequently, 
we define the error operator /C by 

(5.1) Rd^ = -2£_d^u + K.U 

where / = J-'f is the "distorted Fourier transform" from Theorem 14.31 Using the expressions for the direct 
and inverse Fourier transform in that theorem we obtain 



We now return to the variable R as the independent spatial variable instead of r as in the previous section. 
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Integrating by parts with respect to ^ in the second expression we obtain 

(5.2) ICfirj) = ( fimdB. ~ 2^d^]HR, OpiO , ^ v)) - 2 (^1 + f{v) 

where the scalar product is to be interpreted in the principal value sense with / e C(j"((0, oo)). Apriori we 
have 

/C : Co°°((0,cx))) -> C°°((0,cx))) 

therefore we can write 

/>OC 



for a distribution valued function rj —^ K(r],S,)- We refer to (|5.1|) as the transference identity to indicate that 
we are transferring derivatives from R to ^. To asses its usefulness we need to understand the boundedness 
properties of the operator /C. We begin with a description of the kernel Kijj,^). 

Theorem 5.1. The operator JC can be written as 

where the operator /Co has a kernel Ko{ri,S^) of the forrn^'^ 

(5.4) Ko{v,0 = f^Fitv) 

with a symmetric function F(r],(^) of class in (0, oo) x (0,cxd) satisfying the hounds 

1 i + ri<l 



\F{^,V)\< 

\d^F{i,lj)\ + \dr,F{^,V)\< 



sup \did':^F{i,i^)\ < 



where N an arbitrary large integer. 

Proof. We first establish the off-diagonal behavior of K, and later return to the issue of identifying the 
(5-measure that sits on the diagonal. We begin with 1)5. 2|l with / G C^((0,oo)). The integral 



Jo 

behaves like at and is a Schwartz function at infinity. The second factor (f>{R, rf) in (|5.2|l also decays like 
i?2 at but at infinity it is only bounded with bounded derivatives. Then the following integration by parts 
is justified: 

ilKfir^) = (u,C^{R,r^)) = (Cu, cj,{R,^)) 

Moreover, 

Cu - / fiO [C, RdnMR, OpiO + fiO {R9r - 2^d^)mR, OpiO d^ 
Jo Jo 

= / f{0[C,RdB]^{R,OpiOd^+ ^fmRdB-2^d^)cl,{R,Op{Od^-2 Um{R,OpiOd^ 
Jo Jo Jo 

with the commutator 



10. 



The kernel below is interpreted in the principal value sense 
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Thus, 

POO poo 

Lu= f{OW{R)<j,{R,Op{(,)d^+ UmRdR-'2id^)m.i)p{0d(, 
Jo Jo 

Hence we obtain 

7?/C/(77) - /C(C/)(77) = ( r fiOWiRMR, OpiO d^, m V) 







The double integral on the right-hand side is absolutely convergent, therefore we can change the order of 
integration to obtain 

~ OK{v, = p{0(w{R)cf>{R, 0, <P{R, v)] 



This leads to the representation in (|5.4II when ^ ^ rj with 

F{^,r,)^(wiRmR,0,<f>{R,v) 



It remains to study its size and regularity. First, due to our pointwise bound from the previous section. 



(5.5) 



sup|0(i?,e)i 

R>0 


< 








\RdRcl>{R,0\ 


< 


1 3 

min(_R^^ 4 , i?2 


) ve>i 




\d^HR,o\ 


< 


min(i?^^S, i?5 


) Ve>l/2 




|5c</.(i?,0l 


< 


min(i?^ log(l H 


- i?2),r ^|iogeii?) 


VO < e < 1/2 


|5|</>(i?,0l 


< 


mm{R^£,-i,R 






\dlHR,0\ 


< 


min(_R^ log(l H 


- i?2),r ^iiogeii?^) 


VO < ^ < 1/2 



we always have the estimates 

(5.6) \dl^F{^,7j)\<{0-Hvy-' ve + ?7>i 

They are only useful when ^ and rj are very close. To improve on them, we consider two cases: 

Case 1: 1 < ^ + 1], To capture the cancelations when ^ and f] are separated we resort to another integration 
by parts, 

?7F(C,r7) = (w{R)(t>{R,0,C(t>{R,v)) - ([C,W{RmR,0,HR,v)) + CF{^,v) 
Hence, evaluating the commutator, 

(5.7) (?7 - OF{C, V) = -{C^WRdR + Wrr)cP{R, 0, 0(i?, v)) 

Since Wr{0) = it follows that {2WRdR + Wrr)(I){R,(,) has the same behavior as (j){R,0 first region. 

Then we can repeat the argument above to obtain 

iv ' CfFi^, ri) 2WRdR + Wrr]<I>{R, 0,<f>{R, v)) 

The second commutator has the form, with V{R) := —8(1 + R^)^^, 

[C, 2WRdR + Wrr] = AWrrC - AWrrrOr - Wrrrr + ^R^^R-^Wr - Wrr) - 2WrVr - AWrrV 
Since R-^WRiR) - Wrr{R) ^ 0{R^) this leads to 

iv^^fFiS,,^) = (^{W'^^^{R)dR + W^'^^-{R) + ^W^''^'^{RmR,0,<P{R,v)) 
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where by W°'^'^, respectively W^'"'^'^, we have generically denoted odd, respectively even, nonsingular rational 
functions with good decay at infinity. Inductively, one now verifies the identity 

fc-i fc 

j=o e=o 

By means of the pointwise bounds on (j> and dncj) from (|5.5|l we infer from this that 

^fc-j(r;)- | 

Combining this estimate with (|5.6|l yields, for arbitrary N, 

\F{^,v)\<{^ + vyH'^ + \£.^-V^\)~^ provided C + 
as claimed. For the derivatives of F we follow a similar procedure. If ^ and rj are comparable, then from (|5.6(l . 
\97-iF{£^,ri)\ < Otherwise we differentiate with respect to rj in H5.8|l . This yields 

k-l k 

iv - O^'d.Fitv) = ((E^' W^f''{R)dn + J2eW^r{R))<j>{R,^),d,<j>{R,v)) - 2k{7^ - if^"' F{^,j^) 
Using also the bound on F from above we obtain 

<7:— 7T2fe' 1^^'^ 

respectively 
and 

which again yield the desired bounds. Finally, we consider the second order derivatives with respect to ^ and 
rj. For ^ and 77 close we again use the bound from (|5.()|) . Otherwise we differentiate twice in 1)5. 8|) and continue 
as before. We note that it is important here that the decay of W^j"^ and VF^^'^" improves with fc. This is 
because the second order derivative bound at has a sizeable growth at infinity which has to be canceled, 

|9|(/)(i?,0)| « i?5 logi? 

Case 2: ^,?7 ^ 1. Our first observation is that F(0,0) = 0. This is easy to verify by direct integration, 
and is heuristically justified by the fact that W — [L, Rdu]- The pointwise bound 

|5«F(e,77)l <1 

follows by a direct computation. The second order derivative bound is, however, more delicate. We have at 
our disposal the pointwise bounds 

If?7<^<l/2, then these bounds imply that 

\dlr,F{^,v)\< I {R)-^R\\og{l + R^)fdR+ [\ {R)-^ R^ C-^\\og£,\\og{l + R^) dR 

Jo 2 

+ / ,(i?)-2r^'r^|iogci|iog?7|di? 

The main contribution comes from the first term. Integrating we obtain 

|a|,F(C,^)l<|logei' 





3 5 

41] 4 


iv - 






5 




4 


iv- 




e 


3 
4 


iv- 


0'" 
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A similar computation yields, again when rj < ^ < 1/2, 

\dlF{^,r,)\< ( (R)-^RHlog{l + R^)fdR+ (\ log^| log(l + i?^) di? 

Jo 2 
/•OC 

+ / ^(i?)-2^-Jr,3|log^||log77|di?< llogCI' 
It remains to consider the expression dlF{^,T]) for ^ < ry < 1/2. Differentiating in (|^ we obtain 
iv - v) = ^d^F{^, 77) - 0, {2WRdR + Wrr^R, r,)] 



dC 



We differentiate and integrate with respect to rj to obtain 
(5.9) {ri-OdlFitv) 
Using also the bound 



29|.F(e, C) - {dU{R, 0, {2WRdR + WRR.)dccl){R, () 



\dRdcHR,0\ 



< 



i?3 log(l + i?2) R<(-2 

we can evaluate the inner product in (|5.9(l as follows: 



dl^iR,0,{2WRdR + WRR)<l){R,v))\< {R)-'^Rhogil + R^)Rhog{l + R^)dR 



^ (i?)-^i?2log(l + i?2)r^|logC|i?di?+ / ^(i?)-V|logC|i?X~'|logC|i?di?<|logC| = 

5 



Thus, (|5.9() is controlled by 



\{v^Od?,F{^,v)\< / {logCfdC <ii\\ogij\' 



Since ij ^ ^ this yields 

which concludes the analysis of the off-diagonal part of the kernel. 

Next, we extract the S measure that sits on the diagonal of the kernel K from the representation for- 
mula H5.2|) . see also (|5.3|l . To do so, we can restrict ^, 77 to a compact subset of (0, 00). This is convenient, as 
we then have the following asyniptotics of (j){R,S,) for 3> 1; 

Si 



= Re 

{ROr - 2^d^)(j3{R, = -2 Re 



0(i?-") 



3i 



8RC^ 



+ 0{R-^) 



where the O(-) terms depend on the choice of the compact subset. The R~^ terms are integrable so they 
contribute a bounded kernel to the inner product in l|5.2() . The same applies to the contribution of a bounded 
R region. Using the above expansions, we conclude that the J-measure contribution of the inner product 
in H5.2(l can only come from one of the following integrals: 



(5.10) 
(5.11) 
(5.12) 



00 i^OO 



-'0 

2 
1 

2 Jo 



"'O 

00 nOG 



/(Ox(i?)Re 



^a,(«(Or^)«(^)^-^e^^(«'+'") (1 + (1 + 



p{i) didR 



/(Ox(i?) ^95(a(6r')a(^)^"'e-^^(«'-"') ( 1 



'ii 



8Rt] 



8Rr] 



3i 



r 1 + r P(0 d^dR 
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where x is a smooth cutoff function which equals near R = and 1 near R = oo. In all of the above integrals 
we can argue as in the proof of the classical Fourier inversion formula to change the order of integration. 
Integrating by parts in the first integral H5.10|l reveals that it cannot contribute a 5-measure. Discarding the 
R^^ terms from (|5.11|l and (|5.12l) reduces us further to the expressions 



(5.13) 
(5.14) 



OC pOC 



^0 

oo nOO 



3 
8 7o 



/(e)x(i?) Re 
/(Ox(i?) Im 



m^{m~-'Hr])r]--~ 



p(0 d^dR 

R-\r^ -V-^)piOd^dR 



The second integral (|5.14() has both an R^^ and a ~ '7^^) factor so its contribution to K is bounded. 
The first integral H5.13|l contributes both a Hilbert transform type kernel as well as a (5- measure to K. By 
inspection, the 6 contribution is 



Re 



dR 



-n Re 



= -2<2p(^)Re 
-1 , ^p'iO 



-2 piO 

where we used that p{£,)~^ = 7r|ap in the final step. Combining this with the (5-measure in (|5.2|l yields (|5.3|l . □ 
Next we consider the mapping properties for /C. We introduce the weighted spaces " with norm 



l/ll 



(5.15) 
Then we have 

Proposition 5.2. a) The operator /Co from (|5.3|l maps 



-2,Q+l/2 

h) In addition, we have the commutator hound 

Both statements hold for all a Cz R. 

Proof, a) This is equivalent to showing that the kernel 

With the notation of the previous theorem, the kernel on the left-hand side is 



c - 77 

We first separate the diagonal and off-diagonal behavior of Kq, considering several cases. 
Case 1: (^,77) G Q := [0,4] x [0,4]. 

We cover the unit interval with dyadic subintervals Ij — [2^^^, 2-'+^]. We cover the diagonal with the union 
of squares 

2 

j=-oo 
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and divide the kernel Kq into 

^qKq = IahqKq + 1q\aKo 

Case 1(a): Here we show that the diagonal part IahqKq of Kq maps to L^. By orthogonality it suffices 
to restrict ourselves to a single square Ij x Ij. We recall the Tl theorem for Calderon-Zygmund operators, 
see page 293 in \tite\ : suppose the kernel K{rj,£) on defines an operator T : S ^ S' and has the following 
pointwise properties with some 76 (0, 1] and a constant Cq: 

(i) \K{i^,i)\<CQ\^-'n\-^ 

(ii) \K{^,S)^K{t^',0\<Cq\ti-v'\''\^-v\-'-^ for all < 

(iii) \K{T^,i)-K{j^,i')\<CQ\£,~m£,^il\-^-^ forall|e-e'l < 

If in addition T has the restricted boundedness property, i.e., for all r > and ^o,'7o G I^, I|r('^'^''^°)ll2 < 
CqA and ||T*(w'^''"')||2 < CqA where uj''^^°{Cj = lo{{^ - ^o)/r) with a fixed bump-function w, then T and T* 
are L^(IR) bounded with an operator norm that only depends on Cq. 

Within the square Ij x Ij, Theorem 15 . II shows that the kernel of Kq satisfies these properties with 7=1, 
and is thus bounded on L^. 

Case 1(b): Consider now the off-diagonal part Iq\aKo. In this region, by Theorem 15. II 

\KQ{ri,0\< ^ 



/^|log^||log7y| 

which is a Hilbert-Schmidt kernel on Q and thus bounded. 

Case 2: {£,,7]) G Q'^. We cover the diagonal with the union of squares 

00 

B=\Jl,x I, 

and divide the kernel Kq into 

^q-'Kq = IsnQ-'Ko + 1qc\bKo 
Case 2a: Here we consider the estimate on B. As in case la) above, we use Calderon-Zygmund theory. 
Evidently, |it^o(??, 01 ^ ~ on B by Theorem 15.11 To check (ii) and (iii), we differentiate Kq. It will 
suffice to consider the case where the derivative falls on F(^, 77). We distinguish two cases: if — rys | < l^ 
then \£, — rjl < (^"s which implies that 

if, on the other hand, _ jyl] > 1^ then 



which proves property (iii) on B with 7 = 5, and by symmetry also (ii). The restricted property follows 
form the cancelation in the kernel and the previous bounds on the kernel. Hence, Kq is bounded on B. 
Case 2b: Finally, in the exterior region Qc \ B we have the bound, with arbitrarily large N, 

which is bounded by Schur's lemma. 

b) A direct computation shows that the kernel Kq°"^ of the commutator [^9^ , Kq] is given by 

= ivdv + mMv, + KQiv, = ^F-™(C, ry) 
interpreted in the principal value sense and with F"^""^ given by 
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By Theorem 15.11 this satisfies the same pointwise ofF-diagonal bounds as F. Near the diagonal the bounds 
for F'^°"^ and its derivatives are worse^^ than those for F by a factor of (1 + ^) 2 . Then the proof of the 
commutator bound is similar to the argument in part (a). □ 

6. The final equation 

To rewrite the equation (|3.2|l in a final form, we begin by expressing the operator ROr in terms of the 
kernel /C in the transference identity (|5.1|l . We have, with as in Theorem 14.31 



which gives 

This leads to a transport type equation for the Fourier transform x{t, ^) of e by applying J-' to (|3.2|l . Indeed, 
in view of the preceding 



- e2fe-i 



We want to obtain solutions to (|6.1|) which decay as r ^ 00, which means we need to solve the equation 
backward in time, i.e., with zero Cauchy data at r = 00. We treat this problem iteratively, as a small 
perturbation of the linear equation governed by the operator on the left-hand side. For this we need to solve 
the following transport equation 

(6.2) Aidr-2^id^\d\^{T,0^Kr,0, 



A 

We denote by H the backward fundamental solution for the operator 



and by H{t, a) its kernel, 

x{t) 



H{T,a)f{a)da 



The mapping properties of H are described in the following result, which will be proven in the next section. 

Proposition 6.1. For any a > there exists some (large) constant C = C{a) so that the operator H[T,a) 
satisfies the hounds 

c 

(6.3) ||i/(r,(7)||^.,„_^.,„+i/. <rf'^ 

A. 



< 



c 



(6.4) ||(a^__l2^a^)i/(r,a) 
uniformly in a > t. 

This leads us to introduce the spaces with norm 

:=supr^||/(r)||^.,„ 

T>1 

Then the above proposition immediately allows us to draw the following conclusions: 



"'^^The one derivative loss can be avoided by a more careful analysis, but this does not seem necessary here. 
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Corollary 6.2. Given a > 0, let N be large enough. Then 

+ + (dr - 2^^d^)Hb _ < Co ^\\b\ 



with a constant Cq that depends on a but does not depend on N . 

The small factor N^^ is crucial here for our argument to work. On the other hand, the nonlinear operator 
N2k~i from (|6.1() has the following mapping properties: 

Proposition 6.3. Assume that N is large enough and ^ + | > a > ^. Then the map 

is locally Lipschitz from loo,n-2j2,c,+i/2 ^00,^^2,0^ 

The above two results, combined with Proposition 15.21 allow us to use a contraction argument to solve 
equation The next two sections are devoted to proving Propositions 16. II 16.31 Finally, in the last section 

we close the argument. 

7. The transport equation 

Here we consider the backward fundamental solution H for Ht).2|l and prove Proposition 16. II Observe that 
(I6.2|l implies 

[dl + \-\T)i]x{T, \-\t)0 = 6(r, A-2(r)C) 

We introduce the operator 

and the fundamental solutions S'(t, cr, ^), C/(t, cr, ^), which satisfy 

L4,,5(T,a,0 = 0, 5(t,t,0 = 0, drS{T,a,0\r=. = -1 

L5,,C/(t, a, = 0, C/(t, T, C) = 1, drU[T, a, 0\r=a = 

Then H6.2|l may be solved via 



x(r, A-^ (r)C) = - j 5(t, a, 06(fT, A"^ da 
Given this representation, we note that the index a plays no role in H6.3|) and H6.4fl since 



(l + A-2(aK)" 

Hence without loss of generality we set a = 0. Similarly, we can neglect the measure of integration p{^)d£^ 
which also has a polynomial behavior both at and at infinity, 



(7) 



Then the bounds l|6.3|l and l|6.4|l reduce to proving that 

\SiT,a,0\<a(^^y {l + X-\r)0-K |5.5(r, a, ?)| < (^) , l<r<a 

Recalling that A(t) = t^+''' (we are ignoring a multiplicative constant here), we strengthen the first bound 
and prove instead that 

(7.1) \S{T,a,0\<<j[^Y {l + T-i^-K \drS{T,a,0\<{^Y < r < a 

The advantage of doing this is that the last bound is scale invariant. Precisely, one verifies directly the scaling 
relation 

5(T,a,o = c^^(Tr^f^r^i) 
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which leaves H7.1(l unchanged. Hence in what follows it suffices to prove (|7.1|l in the case ^ = 1. We begin by 
constructing two special solutions for the operator Li^r- For small"'^^ r we use a standard WKB ansatz. 

Lemma 7.1. a) (Large r solutions) IJ v is not an even integer then there exist two analytic solutions <j)Q and 
01 of Li^ripj = with a series representation 

oo 

^0 (^) = ^J*^ '^^^ ~ ' = 1 

k=Q 

which is convergent for allr > 0. If v is an even integer then the result still holds with a modification in the 
expression for <j)i , namely 



(/>i(r) =ci(/)o(T)logr + ^CifeT^ , cio = 1 



fc=0 

h) (Small T solutions) There is a solution 02 for Li,t- of the form 

with a smooth and satisfying a(0) = 0. 

Proof, a) We substitute the formal series in the equation 

in the equation and identify the coefficients of the similar terms. This yields 

2k 2k 
Cj,fc(j - - 1 - — ) + Cj.fc-i =0 fc > 1 

Hence the coefficients Cjk can be iteratively computed and satisfy a bound of the type 



< 



(fc!)2 

which implies that the series converges for all r. 

li j = then the argument works for all > 0. If j = 1 then there is an obstruction if v is an even integer; 
indeed, this happens precisely when 2k — v. As usual, this is compensated for by adding in the logarithmic 
term, since 

Li^riM-r) logr) = + T-^2dr(t)n 

has a nonzero coefficient on the r^^ term. 

b) In this case, we use the usual WKB-ansatz which we now recall in a more general setting: we wish to 
solve the equation (9^ + Q)%1} = where Q{t) is a smooth potential for r > 0. Fix some (small) tq > 0. WKB 
means that we seek a solution of the form iIj{t) — '4>q{t)[1 + a(T)] with 

V'o(r) -g-^(T)e*^M, S{t)^ f QH^)dcT 

J To 

Since 

we obtain the following equation for a(r): 
which we solve in the form 



«(t) = - / ij,^{T')i;Ua)Via)[l + a{a)] dadr' 

Q~H<j) [1 - e2*(^^M-^^(-))] V{<j)[l + a{a)] da 



"'^■^The reader should bear in mind that by this r we mean the rescaled one, i.e, ^ 2 t, which can be arbitrarily close to zero. 
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provided these integrals converge at zero. They do in our case: in fact, Q(t) = A~^(r) which imphes that 



a(r) = ci 



r a-^+i [1 - e2-(-" " ')][! + aia)] da 
Jo 



or, after changing variables to a(T'^) = air), 



(7.2) a(T) = iciy 



Jo 



By the boundedness of the kernel, this Volterra equation has a solution a G C([0,oo)) which is clearly then 
also smooth for all r > 0. We now claim that in fact a G C°°([0, oo)). Indeed, the zero order iterate here is a 
smooth function at r = 0: 



for any positive integer m by repeated integration by parts. One now proceeds to show the same for the 
higher Volterra iterates; alternatively, we insert the ansatz 

m 

into (|7.2() and solve for the cocflicicnts dj. In cither case, the conclusion is that H7.2|l has a smooth solution, 
as claimed. □ 

We now use this lemma to prove (|7.1|l . which will then conclude the proof of Proposition 16. II Considering 
the limits at infinity, respectively at 0, one finds that 

This allows us to express the backward fundamental solution S'(t, a) in terms of these bases. Note that we 
suppress the ^ variable as ^ = 1 is fixed. We consider two cases. 
Case 1: a > 1. Then we have 

5'(r, cr) = 0i(cr)0o(T) -0o(cr)0i(T) 

If 1 < T < (7, then (|7.1|1 follows directly from the properties of and ^i. If r < 1 then we express 4'o{t) and 
4>i{t) in terms of the {(/)2, 4>2\ basis to obtain 

^(T,a)= Re(c(a)02(T)), |c(a)| < a 

This gives 

\S{T,a)\ <aT^+^, \drS{T,a)\ <aT-^-^ 

Again H7.1|l follows. 

Case 2: cr < 1. Then we express S{t, a) in the {02, 02} basis to obtain 

^(T,a)= Im (02(^)02 M) 

This gives the bounds 

|5'(t,(t)| <a'^+^T^+^, \drS{T,a)\ <ai+^T-^-^ 

which imply (|7.1|) . 
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8. The nonlinear terms 

In this section we consider the nonlinear source terms, i.e., those given by the right-hand side of l|8.2(l . and 
prove Proposition 1^21 RecaUing that R = rX, we write 

cos{2u2k-i) - cos{2Q) sm{2u2k-i) coa{2eR-i) - 1 
A K^I\2k-i[R = ^ 2e+ — -i 

(8.1) , 

,sm(2eR-2) -2eR-2 
+ cos(2u2fc-i) 



2i?i 

where the regularity of the coefficients above is computed as in Step 2 of the proof of Theorem 12. II 

(8.2) cos(2.,,_0 - cos(2Q) ^ j ^^R-^ ^log Rf , Q._,) 

(8.3) ^^^^^ll^ e/5°(i?- log i?,Q._0 

(8.4) cos(2u2fc-i)e/^°(l,Qfc_i) 

where we used here that tX{t) x r and also that R ^ t (recall the algebras Q and Qk from Definition 12. 
Proposition 16.31 amounts to proving multiplicative estimates in the context of the classical Sobolev spaces. 
Here we use Sobolev spaces adapted to the operator £, namely 

Restating Proposition K . 31 with this notation shows that we need to prove that the map 

e ^ X-^RiN2k-i{R~^e) 

is locally Lipschitz from L<^,N-2jja+i/2 L^,Njja^ ^^^^ ^YiaA has an explicit gain of t'"^ which 

explains why we can improve the time-decay of the first (linear) term on the right-hand side of (|8.1() from 
iV — 2 to A^. On the other hand, there is no such gain in (|8.3|l and H8.4f) . What saves us here is that both the 
second and third terms on the right-hand side of H8.1|) are truly nonlinear terms in e. 
As a technical tool we introduce an inhomogeneous Littlewood-Paley decomposition 



j*oo 

X=l A 



corresponding to a smooth partition of unity {p\} in the Fourier space. Here A G {2-'}°^g and px is adapted 
to frequencies of size A. Our first result is 

Lemma 8.1. Let q e 5(1, Q) and \a\ < ^ + |. Then 

hf\\H^ < ll/llif" 

Proof. We decompose the multiplication operator into its Littlewood-Paley pieces: 

A, /I 

The diagonal sum corresponding to A x is estimated using only the L°° bound on q. For the off-diagonal 
component it suffices to show rapid decay. In fact, we claim that 

\\Px qP^.\\L-^L- < (a* + A)-^-f [log(^ + A)]" 
where m is some large integer. The Fourier kernel of Pa qPf^ is 



KxAv.O = Vpi^fpW) pxiOPt^iv) J q{Rm,R)(t>{v,R)dR 

in the sense that 

^/MHPxqPMi)^ I KxAv,Om)^mdi 
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Therefore, the above bound would follow from the pointwise estimate (recall p(^) x ^ for ^ > 1) 



qiR)^i^,R)^{ri,R)dR 



< 



The function q has a symbol type behavior with respect to R except near R — t, where it has a power type 
singularity (1 — a)"^^ ^ a — R/t, possibly involving also logarithms''^. To separate this singularity from the 
behavior at we use a smooth cutoff to split q into (recall that r is a large parameter) 

q = q<T/2 + q>T/2 

The first term is a symbol of order with respect to R. To proceed, we recall the calculations leading up 
to H5.8|l . The main tool there is the following double commutator identity: if ^ 7^ ry and U is a zero order 
symbol, then 



(e - v)'{u{R)^iR, 0, HR, v)) = ( [[u, c] , c] ^(i?, 0, HR, v) 

(8.5) ^ 

= ( ( - ^Uhb^ + 3R-\Ub.r - R-^Ur) + WrrV + Urrrr + 2URyR + WRRR.dR)(l){R, 0, 0(i?, V) 



where the inner products exist in the principal value sense (recall that V{R) — —8(1 + R^) ^). Iterating this 
identity k times yields 

fc-l k 



3=0 i=0 

nbols of order at most 

For 1 + ^ ^ 1 + ?7 this gives 



where q°'^'^ and qf"^" are symbols of order at most —2k with odd, respectively even, expansions around R = 0. 



\{q<r/2iRmR,0,<f>iR,v))\<{i + v)-' 

for all k which is more than we need. 

The second term <z>t/2 can be thought of as a function of a, 

R 

q>T/2{R) ^ qi{a), a= — 

where qi is supported in [^,2] and has a Q type singularity'^ at a = 1. We divide it into a singular and a 
nonsingular component, 

q>r/2 = q>r/2+q>r/2^ 9>r/2 ■= '?>r/2X[|^_^|^^^^^^pi] , 9>r/2 - «>-/2>^ >(4+^) - i ] ' 

where the x's define a smooth partition of unity relative to the indicated sets. For the singular component 
we bound the integral directly using the pointwise bounds on (j){R, ^) to obtain 

J ql,/,{R)c^{R,OHR,v)dR < |J(l-i?yr)''+^|log(l-i?/r)|™ljl^_^l^^^^^^_.j(0-^(77)-^di? 

< iO-Hvr'^r-'-^^ + r,)-t-i [log(2 + ^ + '?)]" 

< + r?)^5-^ [log(2 + e + V)r 
For the nonsingular component, a fc-fold iteration of (|8.5|l yields 

fe-l k 

(8.6) (e-ry)2'=(<z-/2(i?)0(i?,O,0(i?,'7)) = ( [ E + E ^'Cr(^)] '^(^, 0, ^ ^) 

j=o e=o 

with 

2/c-j-l 2k-e 

<fiR) - E d]i+\i^/,iR), qiyiR) = '^m"(^) d'jiq>r/2iR) 

i=0 i=l 



'■^Strictly speaking, there is a multiplicative constant in a = cR/r, but we ignore it 
'^gi also has a nonsingular part, which by a slight abuse of notation we include in g<T/2 
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where the coefficients are rational functions, smooth for all R> 0, decaying at rates 

The logic behind the numerology here is simple: a factor consumes 2j derivatives, so the remaining Ak 
derivatives need to hit either the symbol q'^^/2(^) °^ weight V (the latter leading to the rational functions). 

We show how to apply these formulas for the case of the even weights, the odd ones being analogous. As 
for the derivatives 

it will suffice to consider two extreme cases: when all derivatives fall on the symbol, or all fall on the cut-off 
function. The contribution by the latter to |((?"*!^^2(-^)'/'(^i "^(^i ^)) I is bounded by (ignoring logs) 

[|K-r|x(^+»,>-5] 

< (C + ,,)-2fer-3-(4fc-2^-2i)^-.-|^^ + ^^_|_|+,^^^,_3 ^^^_3 

as desired. The other cases are checked similarly and we skip them. □ 

This allows us to deal with the coefficients in front of the e terms. As remarked above, the r decay for the 
first term in N2k-i comes from the r^^ factor in the coefficient and from the quadratic (respectively, cubic) 
expressions in i for the remaining terms. It remains to prove the following: 

Proposition 8.2. Let a> \ . Then the maps 

(8.7) e^R-^ (cos(2£ii-3 ) - 1) 

(8.8) £ ^ R~i{sm{2iR-i) - 2iR-i) 

are locally Lipschitz from Hp~^^^^ to H". 

The proof will be split up into the following four lemmas. We first obtain a pointwise bound for frequency 
localized L'^ functions: 

Lemma 8.3. For dyadic X>1 we have 

\Pxf{R)\<Xmm{Ri,X-i}\\f\\L2 

for all f e L'^{m+). 

Proof. Using the inversion formula we write 

/•OO 

Jo 

The pointwise bounds for </), 

|0(i?,e)| <niin{i?i,r^} 

and the Cauchy-Schwarz inequality finish the proof. □ 
We also have estimates for the derivative: 
Lemma 8.4. For dyadic X>1 we have 

\dHPxf{R)\<Xmm{Ri,X-i}\\fh. 



\\dRPxf\\L^<X-^\\f\\L^ 

for all f € L"^ 
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Proof. The first estimate follows from the pointwise bounds on dRcf). For the second bound we can integrate 
by parts (justified by the first bound) to obtain 

M\PxfiR)\\h > {^f, f) > WdufWh + l\\R-'f\\h - cWfWh 

which leads to the desired conclusion. □ 

Next we consider bilinear estimates but with a weight that is singular at 0. This suffices in order to 
estimate the quadratic and the cubic terms in the proposition. The logic behind the Lemma 18.51 is the 
following: dividing by i?2 should amount to a loss of 5* on the Fourier side (since the scaling relation is 
i?f 2 =1). Inspection of the following estimates shows that we do indeed lose a combined | weight in ^ on 
the right-hand side. 

Lemma 8.5. Let a > j. Then 

P-/..ll^..i<ll/ll^..ill.ll^.H 

respectively 

\\R-^f9\\H^<\\f\\^..i\\9\\^^H 
for all /, g so that the right-hand sides are finite. 

Proof. We first use the above pointwise bound to obtain an estimate, 

\\R---PxJPx,g\\L-<min{Xi,X2}\\PxJ\\L4Px29\\L- 

This suffices for both of the above estimates provided that the output is measured at frequency a < 
max{Ai, A2}. Indeed, in that case 

Ai ,A2 cr<max(Ai .X2 ) 

< E K^~'>^2\\PxJh\\Px,9\\2+ E Ar'Al||PAj||2||PA.5l|2 
Ai>A2 Ai<A2 

Ai>A2 

which gives the desired bound since a > j. 

For larger a, however, we need some additional decay. For this we compute using integration by parts 

'R~iPxJPx,9-Pah) - {R-"^PxJPx,g,C''C'''Pah] 
= (^C\R-iPx,fPx,9),C-''P.h] 
To justify the integration by parts we observe that near R — Q we have 

PxJ{R) = RhiR^), g analytic 

Then the bilinear form is given by 

R~iPxJPx29 = Rh{R^), q analytic 
which successively implies that (recall CqR^ = 0) 

C\R--^PxJPx,9) = Rh{R^), q analytic 
We claim that we can estimate the left-hand side here in by 

(8.9) ||/:'^(i?-iPAjPA.5)IU^ <min{Ai,A2}max{Ai,A2}1lPAj|lL2|lPA.5lU^ 

Given the above integration by parts, this implies that 



R--PxJPx,9.P.h 



<min{Ai,A2}max{Ai,A2}'=a-'=|lPAi/|lL2|lPA,.g|lL2|lP./i|U2 
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and further 

||P.(i?-5PAjPA,5)||L^ <min{Ai,A2}max{Ai,A2}'^a-'-||PA,/||LHI^A,5llL= 
thus providing the additional decay for large a. 

It remains to prove (|8.9|) . We assume that Ai < A2 and use different bounds depending on whether R is 

small or large. Assume first that i? < A2 ^ . Then we start from 

oc 00 


Next, we claim that 

(8.11) ||/:^[p-i0(P,C)0(P,ry)]|| 1 <A^i 

If true, then combining (|8.10|l and H8.11|l via Minkowski and Cauchy-Schwarz yields 

||£^Xi?-'^A J(P) Pa.5(^))II < A^-Ai||PA,/||2||PA.g||2 

as desired. To prove H8.11|l . consider first fc = 0. Then 

The higher k cases now follow from Proposition ^31 which allows us to write 

R~^{R,^)(j){R,T]) ^ R^q{R'^,£_R^,'nR^), q analytic 
Then, following our previous discussion concerning applications of wc obtain 

C''{R-^iR,0(b{R,v))^ J2 Rk'v'^qimiR^^R^vR^). g^™ analytic 

e+m<k 

which implies (|8.11|l . 

For large R we use the product rule to write 



/:'=(p-iPA,/PA..g)= w!f^{R)d'j,CP,J-d'^C^Px,g 

2i+2j<2k~e~m 

e,m=o,i 

\W!J^{R)\ < i?-2(fc-'-j)+^+™-| 



Then we have 

WCiR-ip.JP.^g)]] . < J2 xt'~'^~'^\\R-id'j,CP,J-d"^aPx,gh. 

L-^iX^ ,00) ^ — ' 

^ 2i+2j<2k-t-m 

e,m=o,i 

We use Lemmas 18.41 18. 31 to bound the first factor in and the second in L^. This gives 

||/:'=(i?-iPAjPA.5)IL,,-i , < A^-Ai||PA,/|UHI^A..g|U^ 

L^(A2 ,00) 

as desired. □ 

Finally, in order to estimate the higher order terms in the Taylor expansion of the sin and cos functions in 
the proposition we also prove a trilinear estimate: 

Lemma 8.6. Let a > 0. Then 

\\R-'fgh\\H^<\\f\\^^+^\\g\\^^+iJh\\H^ 

for all /, g, h so that the right-hand side is finite. 
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Proof. The pointwise bounds above imply the foUowing estimate, 

(8.12) \\R-'PxJPx,9PxMl-^ < niin{Aj a|}||Pa,/||lHI^a.5||lHI^A3/i||l^ 

which again suffices to estimate the output at frequency a < \ :~ max{Ai, A2, A3}. To see this, we write, 
with the summation variables cr, M, Ai, A2, A3 G {2-'}|^0' 

A>M 

On the right-hand side we distinguish the cases Ai < A2 < A3, Ai < A3 < A2, A3 < Ai < A2. We only treat 
the first case, the other two being similar and easier. Thus, we estimate the right-hand side for fixed M as 
follows: 

II E P^./M{P''P>^JP^29Px.h)f^^< Y {^y^WPx^hWl^i E ^hl\\PU\\L4Px.9\\L^^ 

Ai<A2<A3 " A3>M Ai<A2 

A3>M 

'A3\2", 



A3>M Ai<A2 

<M-^'^\\f\\\^Ag\\\^AhrH^ 



The summation with respect to M is trivial. 

For higher frequency outputs we need some additional decay, 

\\P4Rr'PxJPx,9PxM\L-<^-niri{Xf\]}ma^{\uX2,X3}''a-''\\PxJ\\L^^^^ 

This in turn is a consequence of the estimate 



1 3 



^A JPa2 5^A3 /i) 1 1 < min{ A,^ a; } max{ Ai , A2 , A3 f 1 1 Pa J 1 1 1 1 Pa2 .9 1 1 1 1 ^'a3 /i| I 



which is proved in the same manner as H8.9(l . □ 
These lemmas now imply Proposition 18.21 Indeed, we express the cosine-map in (|8.7() in the form 

i?"5(cos(2£i?"5) - 1) = R-h^q{R^^e^) q entire 
The first factor is bounded by 

\\R-^sHH,<\\e\\\^, 



while for q we use its Taylor series together with Lemma 18.61 which shows that as a multiplication operator 
the factor R~^e^ can be bounded by 

(8.13) \\R-'e'\\H^^H^<\\er . 

tip 

Similarly, we write the sine-map from (|8.8() in the form 

i?ri(sin(2ei?-^) - 2eR-^) = R'^ q{R-^ e^) 
For the first factor we apply Lemma 18.51 twice to estimate 



\R-''eHH^<\\R-^eH ^^x\\e\\ < lle"^ 



if" ^ ^fc c+l t o+l ^ 11-11 



while for the q factor we use again (|8.13|) . 
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9. Proof of the main theorem. 
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Here we summarize how to assemble together the elements of the proof. Fixing > ^ we begin with the 
approximate solution U2k-i given by Theorem 12. II and with the corresponding error e2k-i- The index k is 
chosen sufficiently large, depending on //. Apriori both M2fe-i and 62^-1 are defined only inside the cone 
{r <t}. We can extend them to functions with similar regularity supported in a double cone {r < 2t}. This 
extension is done crudely, without any reference to the equation but insuring the matching on the cone for 
all derivatives which are meaningful. 

With these choices for M2fc-i and e2k-i we seek to solve (13.211 backward in r and find a solution e so that 



(9.1) l|e(r)ll ..i <r^-"', II (^^ + t)^» II ^ ^ ^fc 

Here the exponent a is chosen so that 



1 1/ 

The first bound is solely dictated by estimates for the cubic term in the nonlinearity. The second one is a 
consequence of the regularity of 62^-1; namely, e2fe-i has a singularity of type (1 — aY^^ log™(l — a) on the 
cone. Therefore, if a < 1^/2, then for 62^-1 we have the bound 

\\X-'^Rh2k-i{t{T),X~'R)\\H^ <r-^''+^ 

Using the transference identity we recast H3.2|l for e in the form (|6.1|l with x — Te. By virtue of Proposi- 
tions 16.11 l().3l 15.21 we can solve 1)6.1(1 using the contraction principle with respect to the norm 

||a;|| 2c+i + IK^T - ^)a;||roo,N-ir2,o 

Using again the transference identity and Proposition 15. 21 we return back to e, which has the regularity (|9.1() . 
Now eventually we have to return to the original coordinates (t, r) as well as the function e(i, r). For this we 
define the map 

u{R) TuiR, 9) = e^^R-^^u{R) 

where the right hand side is interpreted as a function in expressed in polar coordinates (i?, 6) . It is easy 
to see that this is an isometry 

T : L2(R+) ^ L^{R^) 
Then for the corresponding Sobolev spaces we have 



Lemma 9.1. For any a > we have 



in the sense that if one side is finite then the other is finite and they have comparable sizes. 

Proof. The spaces iJ^(R+) are defined using fractional powers of the operator C. However, we can also define 
them using fractional powers of the operator Cq since the difference £ ~ £q is bounded in and also in any 
H^. This is easily seen if f3 is an integer, and for noninteger values it follows by interpolation. 
Then the conclusion of the lemma follows from the identity 

ATu — TCqu 

which is valid whenever u £ L'^ and Cqu G L^. □ 

To pass from u{t,R), or alternatively u{t,r), to the co-rotational wave map in terms of the ambient 
coordinates of R"^ D S*^, observe that these coordinates are given hy ((> o T{u), where 4> : R^ — > 5^ C M.^^ is 
given by 

(/){pe^^) — {cos p, ship cos 0, sin p sill 9) 
It is then easily seen that ip o T{u) G i?^"+^(R^), interpreted component-wise. We have now constructed a 
wave map on the cone r < i, < t < which is of class H^~^'^~ on the closure of the cone. To get a solution 
on all of R^+^, extend the solution dtu{tQ, .),u{tQ, .) at time t — tQ to all of R^ within the same smoothness 
and equivariance class. Call the corresponding wave map u{t,r). We claim that this wave map extends to 
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(0, to] X and is of class H^~^'^~ until breakdown at time t = 0. Indeed, by finite propagation speed u{t, r) is 
given by u{t,r) on the light cone r < t,0 < t < to- Furthermore, the u does not develop singularities on the 
interval < i < to, as this could only happen outside the light cone, where energy concentration is precluded 
by the equivariance condition. The fact that singularity formation is tantamount to an energy concentration 
scenario is a consequence of |Sh-Tahj . fTa^ . [Tat] , This concludes the proof of Theorem O 
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